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ABSTRACT 

%^ . Follow-up observations of large numbers of gamma-ray burst (GRB) afterglows, facilitated 

^_' by the Swift satellite, have produced a large sample of spectral energy distributions and light 

curves, from which their basic micro- and macro-physical parameters can in principle be de- 
rived. However, a number of phenomena have been observed that defy explanation by simple 
versions of the standard fireball model, leading to a variety of new models. Polarimetry can 
be a major independent diagnostic of afterglow physics, probing the magnetic field properties 
and internal structure of the GRB jets. In this paper we present the first high quality multi- 
night polarimetric light curve of a Swift GRB afterglow, aimed at providing a well calibrated 
dataset of a typical afterglow to serve as a benchmark system for modelling afterglow po- 
larisation behaviour. In particular, our dataset of the afterglow of GRB 091018 (at redshift 
z = 0.971) comprises optical linear polarimetry (R band, 0.13 - 2.3 days after burst); circu- 
lar polarimetry (R band) and near-infrared linear polarimetry (Ks band). We add to that high 
quality optical and near-infrared broadband light curves and spectral energy distributions as 
well as afterglow spectroscopy. The linear polarisation varies between and 3%, with both 
long and short time scale variability visible. We find an achromatic break in the afterglow light 
curve, which corresponds to features in the polarimetric curve. We find that the data can be 
reproduced by jet break models only if an additional polarised component of unknown nature 
is present in the polarimetric curve. We probe the ordered magnetic field component in the 
afterglow through our deep circular polarimetry, finding Pcirc < 0.15% (2(t), the deepest limit 
yet for a GRB afterglow, suggesting ordered fields are weak, if at all present. Our simulta- 
neous R and Ks band polarimetry shows that dust induced polarisation in the host galaxy is 
likely negligible. 
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1 INTRODUCTION 

Collimated outflows in the fonn of jets are ubiquitous, from AGN 
jets driven by accretion of material by supermassive black holes 
in galactic centres to those associated with stellar sources such 
as X-ray binaries and galactic microquasars. A particularly excit- 
ing view of the fundamental physics of jet sources is offered by 
gamma-ray bursts (GRBs), which form the extreme end of the 
energy- and Lorentz f actor parameter space. Since the discover y 
of afterglows in 1997 (Ivan Paradiis et al.lll997l: ICosta et alJlT997h 
we have a broad picture of the origin and cause of GRBs: through 
a catastrophic event (in the case of long-duration GRBs this is the 
core-collapse of a massive star), a jet of highly relativistic mate- 
rial is ejected. In the standard fireball model (e.g. Meszaros & Rees 
1997), the resulting broadband emission detected from X-ray ener- 
gies through to radio waves - the so-called afterglow - is explained 
by the relativistic ejecta colliding with the surrounding circumburst 
medium. The afterglow radiation we detect is consistent with a 
synchrotron emission (e.g. Meszaros & Rees 1992, 1997; Wijers 
& Galama 1999), characterised by a series of smoothly connected 
power laws, with characteristic break frequencies and fluxes. The 
macroscopic properties of shocks are largely understood, and the 
dynamics of the shock created when the relativistic jet hits the 
circumstellar matter can be written down in terms of the explo- 
sion energy, the density (and density gradient) of the medium into 
which the shock ploughs, and the composition of the shocked ma- 
terial. However, oustanding questions remain on the nature of the 
m/crophysics: how are the relativistic particles that radiate the ob- 
served emission accelerated? Where does the magnetic field in the 
shocked region come from and what is its structure? 

^Since its launch in 2004, the Swift satellite (' Gehrels et al.l 

|2004) has provided us with hundreds of well sampled X-ray and 
UV/optical afterglow light curves from seconds to mont hs after 
burst (see for a review iGehrels. Ramirez-Ruiz & Foxll2009h . These 
light curves revea led complexity tha t was unexpected from pre- 
Swift light curves dPiran & Fanll2007h : the standard fireball model 
has received a growing number of modifications, and the concept 
of a canonical light curve was introduced to explain the steep fades, 
plateau phases and multiple breaks in Swift X-ray light curves 
(jNpusek et al. 2006). In particular, discrepancies between X-ray 
and optical light curves of GRB afterglows and the rich array of 
light curve morphology including extended plateaux, rebrighten- 
ings and flares ch allenge existing GRB models (see for a review 
IPiran & Fanll2007h . Explanations for these features include com- 
plex jet structure, energy re-injection due to late-time central en- 
ergy activity, variable microphysics and ofl'-axis emission, with a 
combination of efl'ects likely in action in most afterglows. How- 
ever, these model variations are in principle non-degenerate: linear 
and circular polarimetry have the ability to separate out the various 
models for the behaviour of the (early) afterglows in a completely 
independent way, as these models have difl'erent implications for 
the magnitude and orientation of the optical pola risation, as well as 
for their time-dependence (e.g. lRossi et al.ll2004l) . 

The ground-breaking discovery and interpretation of sudden 
achromatic steepening in light curves at ~ 1 day post-burst in pre- 
Swift GRBs - so-called jet breaks - provide d the first convincing 
evidence of h ighly collimated ejection (eg. iRhoadsl 1 19971 . 1 19991 : 
ISari et al.lll99 9) and allowed jet ope ning angles and collimation- 
corrected energies to be derived (e.g. iFrail et al.luOOlb . Jet breaks 
were expected to be ubiquitous in Swift light curves, but the added 
complexity and multiple breaks in these light curve s have made 
jet breaks difficult to identify unambiguously (e.g. ICurran et al.l 



120071) . In contrast, the linear polarisation around the time of the 
jet break is predicted to have a unique signature (IS aril 119991 
bhisellini & Laz zati"l999': lRossi et"ani200 4'): At early times a dis- 
tant observer located slightly off'-axis will see the afterglow as a 
ring shaped source, which has strong polarisation in radial direc- 
tion (assuming the magnetic field is compressed in the plane nor- 
mal to the motion, i.e. ordered in the plane of the shock). In the 
received light, integrated over the ring, the polarisation will largely 
cancel out. However, as the fireball decelerates the size of the ring 
increases and at some point the edge of the jet is reached, at which 
point the polarisation does not completely cancel out anymore and 
linear polarisation is observed. As the ring expands further, the op- 
posite edge of the jet is reached as well, giving rise to a second 
peak in the linear polarisation curve. Therefore polarimetry pro- 
vides a useful tool to unambiguously identify jet breaks from other 
light curve breaks, though measured polarisation curves of GRB 
afterglows have in some cases shown puzzling features (see e.g. 
the cases of GRBs 030329 and 020405 discussed further in Section 
4.1). In addition, the linear polarisation properties probe jet inter- 
nal structure, as different polarisation properties are predicted for 
'uniform' or 'structured' jet energy distributions (Rossi et al. 2004, 
see also Section 4.2.2). 

It has been shown that if some fraction of the magnetic field 
in the shock has large-scale order, the light may also be circu- 
larly polarised (Masumiya & loka 2003; Sagiv, Waxman & Loeb 
2004). For the forward shock (in the ambient medium), most mod- 
els predict 0.01-1% circ ular polarisation during the first day (e.g. 
iMatsumiva & Iokall2003h . depending on the strength and order of 
the magnetic fields and the wavelength of observation. Part of this 
range is within reach of large telescopes, as we will demonstrate in 
Section 3.2. 

In the Swift era few polarimetric measurements have been 
performed. Arguably the most successful have been the early 
time linear polarisation measurements of the afterglows of GRBs 
060418 and 090102, performed with the RINGO polarimeter on 
the rapidly responding robotic Liverpool Telescope (iMundell et al.l 
I2OO7I : ISteele et al.l l2009), the latter GRB s howing a linear po lari- 
sation of 10% just 190 seconds after burst. ISteele et all (l2009h in- 
terpret the light curve properties and high degree of polarisation 
in GRB 090102 - which was detected when emission from the re- 
verse shock dominated the received optical light - as evidence for 
large-scale ordered magnetic fields in the expanding fireball. This 
detection of significant polarisation at early times provides strong 
motivation for following the polarisation properties of the fading 
afterglow for many hours or days after the GRB. Howeve r, late- 
time polarimetric observations of GRBs (see lLazzatill2006l for an 
overview) show that the combination of fading afterglow and low 
levels of polarisation (few percent) require large telescopes to ob- 
tain polarimetric observations over multiple days to high precision 
(<0.3%), using multiple filters to distinguish dust-induced from in- 
trinsic afterglow polarisation. In this paper we present the first such 
dataset for a Swift burst, and the most extensive polarimetric dataset 
of any burst since GRB 030329 (Greiner et al. 2003). 

This paper is the first in a series on the afterglow of 
GRB 091018. In this first paper we will describe our dataset in de- 
tail, and compare the observed polarisation light curves with those 
of other {^vq- Swift) bursts and theoretical models. In a forthcoming 
second paper (Wiersema et al. 2012 in prep; hereafter called Pa- 
per 2) we will discuss the dust properties in GRB sightlines from 
a combination of multi- wavelength polarimetry, spectroscopy, host 
imaging and broadband afterglow spectral energy distributions. In 
a third paper (Paper 3) we will model the afterglow physics of 
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Figure 1. Shown are small sections of dataframes from the ISAAC Ks band 
polarimetry (Left; angle frame from linKl, see Section 2.2) and optical 
FORS2 polarimetry (Right; angle -45 from circular polarimetry epoch cirl). 
The circles are 5 arcseconds in radius in each image, and show the source 
in the o and e beams. 



the brightness of the afterglow we switched to circular polarimetry, 
for which the models predict a much lower degree of polarisation, 
requiring a large number of detected photons. Directly following 
the circular polarisation, we continued our linear polarimetry. We 
obtained a further 6 sequences, and a further 2 at the end of the 
night. The next night we obtained further datasets at the beginning 
and the end of the night. A last, deep, point was obtained in the 
third night. Note that for all F0RS2 polarimetry we position the 
afterglow close to the centre of the F0RS2 field of view, on chip 1, 
where instr umental polarisation is ver y low (well below our statis- 
tical errors. IPatat & RomanieUoll2006l) . 



GRB 091018 in greater detail and fit the polarisation data to more 
detailed models, in particular considering energy injection. 

This paper is organised as follows: in Section |2] we describe 
the polarimetry, spectroscopy and broadband photometry observa- 
tions, the data reduction and data calibration techniques; in Section 
3 we discuss the detected features in the data; and in Section 4 
we compare the data to models and previous studies of afterglow 
polarimetric light curves. Throughout this paper we adopt a cos- 
mology with Ho = 71 kms-^ Mpc'^ Q,„ = 0.27, ^a = 0.73. 



2 OBSERVATIONS 

GRB 091018 triggered the Burst Alert Telescope (BAT) onboard 
the Swift satellite at 20:48:19 UT on October 18, 2009 (trigger 
373172; Istamatikos et al.ll2009h . The prompt e mission shows the 
burst to likely belong to the class of long bursts (Kouveliotou et al. 
1993h. with a durati on Tgo = 4.4 + 0.6 s ( Markwardt et al. 2009; 
Ukwatta et aDl2012l report a spectral lag of 143 + 297 msec). An 



X-ray and optical afterglow were found by th e Swift X-ray tele - 
scope (XRT) and UV-optical telescope fUVOT). [ch"en et all OOO^ 
reported an afterglow redshift of 0.971 shortly after. Based on the 
initial brightness of the UVOT afterglow we activated our VLT po- 
larimetry programme (programme 084.D-0949, PI Wiersema). 



2.1 R band linear- and circular polarimetry 

2.1.1 Data acquisition 

We acquired imaging polarimetry using the Focal Reducer and low 
dispersion Spectrograph (F0RS2) mounted on Unit Telescope 1 
(Antu) of the VLT, using the F0RS2 Rspeciai filter. After passing 
through a half- or quarter wavelength plate (in the case of linear 
and circular polarimetry, respectively), a Wollaston prism splits the 
incoming light into a ordinary and extraordinary beam (hereafter 
the o and e beam) that are perpendicularly polarised. These o and 
e beams are imaged simultaneously, using a mask to avoid over- 
lap of the two beams on the chip, as shown in Figure [T] For each 
linear polarisation data set we used four rotation angles of the half 
wavelength plate, of 0, 22.5, 45 and 67.5 degrees. The circular po- 
larisation data were taken with -45 and +45 degree angles of the 
quarter wavelength plate. 

To ensure relatively homogenous polarimetric uncertainties as 
a function of time, exposure time was increased as the source faded. 
We employed a small dithering pattern, while taking care not to po- 
sition the afterglow too close to the edges of the mask. We acquired 
a total of 20 linear polarisation series and 4 circular polarisation 
series, see Table [T] 

We began our F0RS2 polarimetric monitoring with a single 
short sequence of linear polarimetry. After this, to make full use of 



2.1.2 Data reduction and calibration 

We bias and flat field corrected all data using tasks in IRAF and 
using standard twilight sky flats. We use our own software in com- 
bination with IRAF routines to perform aperture photometry on the 
o and e images of the afterglo w to find the sou rce fluxes fo and 
/e, using the same approach as IRoI et alj (l2003h . We used a see- 
ing matched aperture of 1.5 times the on-frame full width at half 
maximum (FWHM) of the point spread function, fitted per frame 
independently for the o and e as small diff'erences in PSF shape 
may occur between o and e beams, particularly for objects far off"- 
axis. The sky subtraction was done using an annulus of inner and 
outer radii 3 and 4 times the FWHM. We measure fo and fe for all 
point-like objects in all frames. 

We express the polarimetry information in terms of the Stokes 
vector T = (2, t/, V,/) (see e.g. Chandrasekhar 1950), where the 
components of this vector have the following meaning: Q and U 
contain the behaviour of the linear polarisation; V the circular po- 
larisation and / is the intensity. Generally we will use the nor- 
malised Stokes parameters q - Q/I, u = U/I and v = V/I in 
this paper. Theoretical models are often expressed in terms of the 
polarisation degree P and polarisation angle 6. These quantities can 
be found from the Stokes parameters through the relations: 



P]m — 



^[^TU^ 



6= -arctan(t//2) + 

Par = V/I 

Constant is an off'set defined so that the resulting angle 6 conforms 
to the standard definitions of position angle (angle from North, 
counterclockwise, see Fig|2]): 

r 0° if2>0andt/^0; 
= i 180° if2>0andt/<0; 
[ 90° if 2 < 0. 

Note that the conversion from Stokes parameters to P brings 
with it complications, discussed further below, so wherever possi- 
ble we will work in Stokes parameter space. 

To measure the Stokes Q/I and U/I parameters a measure- 
ment at two rotation angles suffices, but this introduces substantial 
systematic uncertainties, caused by imperfect flatfielding, imper- 
fect background subtraction and imperfect behaviour of the half 
wavelength plate and Wollaston prism. Using four rotation angles 
with a constant stepsize (in this case of 22.5 degrees) means that 
several of these systematic eff'ects cancel out (in particular the ef- 
fects of background subtraction and flatfielding), leading to greatly 
improved polarimetric error terms. With these angles one can write 
(Patat & Romaniello 2006): 
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Table 1. Log of our polarimetric observations. The ID column gives a label to each polarimetric dataset to make discussion easier, with linear polarimetry 
datasets labelled "lin" and circular polarimetry "cir". The exposure time is the exposure time per retarder angle (4 angles for linear, 2 for circular polarimetry). 
The QII and f/// values are as measured: to obtain values corrected for scattering on Galactic dust, subtract the values derived in Section 2.1.3. 
*: combined value of three complete series (see Section 2.2). 



ID 



FORS2, R. 



-special 



T-To 
(mid; d) 



Exp. time 

(s) 



Q/i 



u/i 



v/i 



linl 


0.13189 


180. 


cirl 


0.14302 


300. 


cir2 


0.15133 


300. 


cir3 


0.15964 


300. 


cir4 


0.16795 


300. 


lin2 


0.18043 


240. 


linS 


0.19701 


300. 


lin4 


0.21282 


300. 


lin5 


0.23030 


300. 


lin6 


0.24611 


300. 


lin? 


0.26357 


300. 


linS 


0.27936 


300. 


lin9 


0.45483 


300. 


linlO 


0.47064 


300. 


linll 


1.13940 


300. 


linl2 


1.15522 


300. 


linl3 


1.17345 


300. 


linl4 


1.18926 


300. 


linl5 


1.20729 


300. 


linl6 


1.22310 


300. 


linl 7 


1.36006 


600. 


linlS 


1.39183 


600. 


linl9 


1.44926 


600. 


linlO 


2.39023 


720. 



-0.0014 ± 0.0009 



-0.0036 ± 0.0009 



-0.0005 ± 0.0014 
0.0011 ±0.0015 
0.0003 ± 0.0015 

-0.0017 ± 0.0015 



0.0009 

-0.0012 

0.0010 

-0.0006 

0.0083 

0.0056 

0.0057 

0.0119 

0.0067 

-0.0162 

-0.0167 

0.0088 

-0.0123 

-0.0070 

-0.0126 

0.0022 

0.0014 

-0.0141 

0.0111 



± 0.0009 
+ 0.0008 
± 0.0009 
± 0.0009 
± 0.0009 
± 0.0009 
± 0.0009 
± 0.0012 
± 0.0012 
+ 0.0023 
± 0.0022 
+ 0.0023 
± 0.0023 
± 0.0024 
± 0.0024 
±0.0018 
±0.0018 
± 0.0019 
± 0.0023 



-0.0040 

0.0024 

-0.0023 

-0.0037 

-0.0039 

-0.0025 

-0.0063 

-0.0025 

-0.0009 

0.0078 

0.0260 

0.0130 

-0.0147 

-0.0054 

-0.0080 

-0.0056 

0.0057 

-0.0031 

-0.0092 



0.0009 
0.0008 
0.0009 
0.0009 
0.0009 
0.0009 
0.0009 
0.0012 
0.0012 
0.0023 
0.0023 
0.0023 
0.0023 
0.0024 
0.0024 
0.0018 
0.0019 
0.0019 
0.0023 



ISAAC, Ks 



linK* 0.4309 



720.* 



0.0204 ± 0.0078 



0.0008 ± 0.0083 



2 2 

Q - QII - ITrZi^i ^^^ (i^/^) ^rid u = U/I - t; y,P^i sin (in/l) , 



NORTH 



where N is the number of half wavelength plate positions, and Ft is 
the normalised flux diff'erence of a source in the o and e beams at 
the i-ih angle 

Fi = (fo,i - fed I (fo,i + /..) = (fo,i - fed //. 

For the circular polarimetry, taken with two angles of the retarder, 
we can simply write 

v = y//=i(F45-F_45). 

After measurement of the fluxes of the sources in the images 
we compute their observed Stokes parameters and their uncertain- 
ties using the expressions above. The afterglow is always posi- 
tioned on nearly the same position on CHIPl (barring a small dither 
of ~ 16 pixels in Y direction in some epochs), on the optical axis. 
As such we expect the instrumental linear polarisation to be negli- 
gible; see Patat & Romaniello (2006) for a thorough error analysis 
of FORSl, whose polarimetry optics were moved to the F0RS2 in- 
strument in 2009. Additionally, circular polarimetry Stokes V/I has 
no instrumental polarisation on axis. In conclusion, the values of 
the Stokes parameters found through the methods above allow an 
investigation of the time dependence of the afterglow q, u and v. 

Theoretical models are generally expressed in terms of the lin- 
ear and circular polarisation fractions Pun and Pcir, and we therefore 



EAST— 4 
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Figure 2. Coordinate definitions used in this paper. 



convert the Stokes parameter information to these quantities us- 
ing the equations above. We note that the uncertainty on the linear 
polarisation angle is a function of the intrinsic polarisation degree 
(cTe = CTpj.^ /2Piin), which means that for the low polarisation values 
and faint fluxes of afterglows the uncertainties in 6 are very large. 
Finally, the position angles found are corrected for the instrumental 
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zero angle offset, which is -1.48° at 655 nm (the central wavelength 
ofthefilterfl 

Errors on q and u are generally distributed as a Normal dis- 
tribution and the Stokes parameters can have positive and negative 
values. I n contrast, Pun is a positive definite quantity. As demon- 
strated in lSimmons & Stewar t* ('1985') (their Section 2), integrating 
the equation for the distribution of measured (P, 6) over 6 shows that 
Piin follows a Rice probability distribution rather than a Gaussian 
one. Directly using the equations above will lead to overestimated 
Piin and incorrect confidence intervals, an eff'ect generally referred 
to as the linear polarisation bias. The correction to the resulting Pun 
and associated confidence ranges has been studied through both an- 
alytical and numerical (Monte Carlo) methods. Generally speaking, 
this correction depends on crp, P, and the signal to noise (SNR) 
of fe + fo (i.e. the SNR of It). In the literature one often sees the 
Wardle & Kronberg (1974) prescription used, in which the input P 
values are multiplied by ^1 - (crp/Py to find the bias-corrected 
polarisation. We follow Sparks & Axon (1999) in using a parame- 
ter 77 = P x SNRi7 to trace the expected behaviour of the bias and 
(Tp: when 77 > 2 the Wardle & Kronberg correction is valid, and crp 
is as computed directly from the uncertainty of F. In our F0RS2 
dataset before correction for Galactic polarisation (Section [2. 1.31) 
all data points have 77 > 2, and the bias correction is in all cases 
small compared to crp. 



2.1.3 Polarisation induced by the Galactic interstellar medium 

To bring the afterglow q and u to an absolute frame, we need to 
correct for the linear polarisation induced by scattering on dust in 
our Galaxy (Galactic interstellar polarisation, GIP). Note that VII 
does not require a correction as there is no circular GIP and the 
resulting linear polarisation is much smalle r than the threshold for 
noticeable linear - circular crosstalk (Patat & Romanielloll2006n . 

We assume here that the average intrinsic polarisation of a suf- 
ficiently large sample of field stars is zero, that the bulk of the po- 
larising dust is between these stars and the observer, and that there- 
fore the observed distribution of q and u of field sources measures 
the GIP. We measure the Stokes q and u values of all sources that 
have a SNR above 20 in the F0RS2 field, on both chips, and in all 
epochs. A large number of artefacts are visible in our F0RS2 imag- 
ing polarimetry, which have the appearance of a group of stars, and 
which are likely caused by reflections at the retarder platej^; we take 
special care to avoid stars close to these artefacts. 

We then eliminate sources that appear significantly extended 
(galaxies), which may give spurious linear polarimetric signals if 
seeing conditions change during a polarimetric sequence. We fur- 
ther eliminate sources that are close to the mask edges, and sources 
that are close to saturation in one or more of the images. 

It is well known that FORSl, whose polarimetric optics are 
now cannibalised in F0RS2, displayed a roughly radial instrumen- 
tal linear polarisation pattern, with polarisation magnitude depend- 
ing on the radial distance to the optical axis (see Patat & Ro- 
maniello 2006). The shape of the linear polarisation pattern and 
its magnitude depend on the wavelength, and was found to reach 
nearly ~ 1.5% at the edges of the field of view (Patat & Romaniello 
2006). While this eff^ect is well calibrated for FORSl in V and / 
band, a similar analysis has not yet been performed for F0RS2. We 
therefore use a conservative approach and choose only stars close 



(radially) to the GRB position. It is of some importance to not select 
only the brightest stars (which would give smallest uncertainty in 
the GIP value), as this may result in a bias towards the nearest stars 
(in distance), and therefore may not sample the complete Galactic 
dust column towards the GRB. We also make sure stars are selected 
with a fairly uniform distribution over the field (in terms of position 
angle). 

We determine the GIP q and u values using three methods: 
we find the weighted mean of the field star distribution; we fit 1- 
dimensional Gaussian curves on the histograms of the q and u dis- 
tributions; and we fit a 2-dimensional Gaussian function (with its 
orientation as a free parameter) on the {q, u) distribution after con- 
verting this to a 2-dimensional histogram (where a range of bin 
sizes are used based on the number of input data points). We per- 
form these fits for a range of limit values on the maximum radial 
distance r between the input stars and the GRB; and on the maxi- 
mum allowed crp of the stars. Unfortunately there are only a small 
number of fairly bright stars near to the GRB position, making it 
necessary to use a large limit on the radial distance to obtain a suf- 
ficient number of datapoints. We find a best balance between the 
number of sources in the fit and the uncertainties on the resulting 
^Gip and wgip for limits crp < 0.35% and r < 2.0 arcmin (103 dat- 
apoints). Note that r - 2 arcmin would correspond to ~ 0.23 and 
~ 0.34% i nstrumental polarisat i on in V and / band in FORSl, re- 
spectively; |Patat^Roinaniello| (120061) . Using these limits the three 
methods give very similar GIP values. In Figure |3] the distribution 
of datapoints (103 points) is shown. We find ^gip = -0.0028 and 
uqiv - -0.0036, or Pgip = 0.45%, ^gip = 116°. The 2-dimensional 
fit gives the width of the Gaussian along the q axis as cTq = 0.0030 
and along the u axis as cr^ = 0.0025, and the rotation of the con- 
tainment ellipse from the q axis in radians, counter-clockwise, as 
= 1.7 radians. The evidence for dev iation from radial symmetry 
is not very strong, but the analysis of IPatat & Romaniello! (120061) 
shows that the FORSl off'-axis polarisation pattern exhibited non- 
radial behaviour in the B band. We use this elli pse as a conserv ative 
error in the GIP Stokes parameters, following R ol et al.l(l2003h . The 
fairly large uncertainty in these is likely dominated by the contribu- 
tion of the unknown instrumental polarisation pattern and the low 
number of bright stars within ~ 1 arcminute of the GRB position. 

The empirical relation Pg ip < 0.09E(B - V) 
dSerkowski. Matheson & Fordlll975h gives Pgip ~ 0.3%, which 
compares well with the value derived above. 

We now subtract from the GRB q and u values the ^gip and 
wgip found above, to remove the (constant in time) component to the 
polarisation caused by the Galactic interstellar medium. We prop- 
agate the errors on the GIP Stokes parameters into the resulting 
corrected values ^GiPcorr and WGiPcorr- As is clear from Table [T] the 
afterglow polarisation at very early times (e.g. linl, linS) all but 
vanishes, indicating that at early times the detected polarisation is 
dominated by GIP. In many epochs, the ^GiPcorr and WGiPcorr errors 
are dominated by the uncertainties in the GIP Stokes parameters. 
We compute GIP corrected values for the polarisation and position 
angle, which will bring the measured values onto an absolute scale. 
It is clear that the errors crp are increased dramatically, while P de- 
creases in several epochs (to near zero in e.g. linl), i.e. a drop in 
P/o-p. At the same time, the statistical detection signal to noise (i.e. 
SNR of //, simply referred to as SNR in Section 12.1.2) stays the 
same. As the polarisation bias is a function of both SNR and P/crp, 
the polarisation bias correction is much larger in several epochs af- 
ter GIP correction, showing the necessity of treating GIP eff'ects 
^ Tabulated at the FORS instrument webpages in q, u space rather than P, 6 space in low polarisation regimes. In 

^ |http : //www ■ eso . org/observing/df o/quality/FORS2/qc/problems/prs )fe^fflia^.tfa>3:e5gtteiWardle & Kronberg prescription is not valid (see 
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Figure 3. Measurements of q and u of field stars to determine the contri- 
bution to the received polarisation by scattering on dust within our Galaxy. 
Dashed lines indicate ^ = and m = 0. Binsize used in this plot is 0.0025. 



[Simmons & Stewart 1985;, and we use the methods recommended 
by [Simmons & Stewart ( 1985) for low Pjcrp situations (which uses 
a maximum likelihood estimator, this estimator has lowest bias fac- 
tor and lowest risk function values in this signal to noise range). In 
two cases this led to robus t zero polarisation values. We refer to 
I Simmon s & Stewarg (Il985h for a discussion on the probability of 
correctly inferring zero polarisation for different estimators. 
All resulting values are listed in Table |2] 



2.2 K band linear polarimetry 

2. 2. 1 Data acquisition 

We obtained an imaging polarimetry data set using the short wave 
polarisation (SWP) mode of the ISAAC instrument mounted on 
Unit Telescope 3 (Melipal) of the VLT, using the Ks filter. Obser- 
vations coincided in part with R band linear polarimetry, see Table 
[T] ISAAC SWP mode uses a similar approach to F0RS2 to obtain 
polarimetry: a Wollaston prism is inserted in the light path, and 
the incoming light is split into a ordinary and extraordinary beam 
(the o and e beam) that are perpendicularly polarised. The o and e 
beams are imaged simultaneously, using a mask to avoid overlap. 
The angular separation of the o and e images is wavelength depen- 
dent. The main difference between ISAAC and F0RS2 imaging 
polarimetry is that a half wavelength plate is not present in ISAAC, 
requiring rotation of the instrument to obtain observations at the 
necessary angles; and the location of ISAAC in the Nasmyth focus. 
At each instrument angle we use integration times of 4 X 15 
seconds, executed 4 times with significant dithering to facilitate sky 
subtraction. Because of the aperture mask, less than half of the field 
of view is visible at any given rotation angle, and one generally 
uses a dither pattern designed to sample the entire field. However, 
as the afterglow is faint and polarisation low, we dither such that 
the afterglow always stays within the mask opening. After these 
4x60 seconds exposure, the instrument is rotated, and observations 
are repeated for the new position angle. By mistake we employed 
cumulative rotation offsets of 0, 22.5, 22.5 and 22.5 degrees. As 
ISAAC does not have a halfwavelength plate, the better choice of 
cumulative angles would have been 0, 45, 45, 45, so as to achieve 



ID 


Plin 


e 


Ain,GIPcon- 


^GlPcorr 




(%) 


(degrees) 


(%) 


(degrees) 


linl 


0.36 + 0.13 


126.1 + 19.3 


(0.32)''' 




linl 


0.38 + 0.12 


143.1 + 17.7 


0.21+0.31 


177.0 + 47.5 


linS 


0.24 + 0.11 


60.2 + 25.1 


0.56 + 0.27 


37.7 + 24.7 


lin4 


0.22 + 0.12 


148.5 + 28.4 


0.26+0.31 


9.2 + 43.7 


linS 


0.35+0.12 


132.1 + 19.0 


(0.32)t 




lin6 


0.90 + 0.12 


168.9 + 7.9 


1.07+0.30 


179.2 + 16.1 


lin7 


0.59 + 0.12 


169.5 + 12.0 


0.78+0.31 


3.9 + 21.2 


linS 


0.84 + 0.13 


157.4 + 8.8 


0.84 + 0.30 


171.1 + 19.9 


lin9 


1.20 + 0.17 


175.5 + 8.1 


1.44 + 0.32 


2.2 + 12.6 


linlO 


0.65+0.17 


177.7 + 14.7 


0.94 + 0.32 


8.0+18.6 


linll 


1.77 + 0.32 


78.6 + 10.2 


1.73+0.36 


69.8 + 11.7 


linl2 


3.07 + 0.32 


62.8 + 5.9 


3.25+0.35 


57.6 + 6.1 


linl3 


1.53 + 0.32 


29.5 + 11.9 


1.99+0.35 


27.6 + 10.0 


linl4 


1.89 + 0.32 


116.4 + 9.8 


1.42+0.36 


114.6 + 14.0 


linl5 


0.81+0.33 


110.2 + 21.7 


0.27 + 0.38 


101.6 + 47.1 


linl6 


1.45+0.33 


107.7 + 13.0 


1.00 + 0.38 


102.2 + 20.2 


linl 7 


0.54 + 0.25 


147.0 + 24.0 


0.41 +0.34 


168.8 + 36.6 


linl8 


0.52 + 0.26 


39.5 + 25.3 


0.97 + 0.32 


32.8 + 17.8 


linl9 


1.42 + 0.27 


97.6 + 10.7 


1.08+0.35 


88.7 + 17.9 


linlO 


1.40 + 0.32 


161.6 + 13.1 


1.45+0.37 


169.0 + 14.3 


linK* 


2.0 + 0.8 


10.9 + 20.9 


2.0 + 0.8 


10.9 + 20.9 



Table 2. Linear polarisation values, where coordinate definitions are used as 
in Figure 12] All values have been corrected for polarisation bias (after cor- 
rection for GIF in the case of the last two columns), see Section [2X3] The 
uncertainties in the GIF q, u value are fully propagated, and dominate the 
uncertainties in the GRB afterglow GIF-corrected q, u. The GIF corrected 
polarisation is very low (with inferred values of zero in two cases), and as 
such the errors on 6 are very large (as ere is a function of 1/P, see Section 

[2X3. 

*: Combined value from linKl, linKl and linKS. The GIF correction in K 

band is more than an order of magnitude smaller than the statistical error in 

the polarisation, so no correction is made. 

^: See Simmons & Stewart (1985) for the probability of correctly inferring 

true polarisation po = from measured polarisations P using a maximum 

likelihood estimator. The number in brackets is the 67% confidence interval 

boundary. 



the same effects as in the FORS data (switching of the o and e 
images to minimize flatfield errors and Wollaston throughput un- 
certainties). These series of observations are repeated three times, 
resulting in complete sets linKl, linK2 and linK3. Note that be- 
cause of the rotation of the instrument, field stars are only visible in 
some of the rotation angles, and rotate behind the mask at other an- 
gles. We note that a series of data were taken before linKl with an 
erroneous dithering pattern: the afterglow disappeared behind the 
aperture mask at the third and fourth angle. In further analysis we 
will only use the complete cycles linKl - linK3: the calibration as 
described below requires measurement of both q and u for a given 
parallactic angle. 



2.2.2 Data reduction and calibration 

Data were corrected for dark current, flat fielded and background 
subtracted using tasks in IRAF, and the 4 exposures of 60 seconds 
per angle per cycle were co- added. Aperture photometry was per- 
formed on the e and o images of the afterglow, in the same way 
as described above for F0RS2. We compute the normalised Stokes 
parameters: 

(h,0 - h,o) J rr/T (4,45 -^e,45) 



Q/i- 



and 



U/I: 



(h,0 + 4,o) (4,45 + 4,45) 

using the data with rotation angle and 45 degrees. As we did not 
take data with 90 and 135 degree angles, we can not use these to 
cancel out systematics as we did with the F0RS2 data. We expect 
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Figure 4. Optical light curves used in this paper (for their values, see Ta- 
ble 4). All magnitudes are in the AB photometric system, and the FORS2 
Rspeciai magnitudes are converted to r. Light curves in different bands are 
offset for clarity. 



these effects to be of the same order or smaller as the statistical 
errors (the afterglow was faint in K). We also compute Stokes pa- 
rameters from the 22.5 and 67.5 degree data: 



q2 = (Q/I)2 



(4,22.5 - ^e,22.5) -, ,tj n\ ih ,61 .5 - h ,61 .s) 

— and U2 = (t///)2 = — r, 

Uo,22.5 + ^e,22.5J Uo,67.5 + h,61.5) 



which we will calibrate independently (they differ from {q, u) by a 
rotation in coordinate basis). 

ISAAC is a Nasmyth focus instrument, and as such less well 
suited to accurate polarimetry: the instrumental polarisation of Nas- 
myth instruments is dominated by reflection of the tertiary mirror 
M3 (at 45 degrees, therefore introducing strong polarisation which 
is highly dependent on position or parallactic angle), which re- 
quires particular care in calibration. Our calibration strategy closely 
follows the recommendations of Joos et al. (2008) and Witzel et al. 
(2011). These authors highlight in their papers the great benefits of 
using a train of Mueller matrices to create a polarimetric model of 
the instrument plus telescope for Nasmyth mounted instruments: 
all elements contributing to instrumental polarisation effects are 
described in 4 x 4 matrices acting on the Stokes vector S . These 
matrices are then multiplied to give a Mueller matrix of the whole 
telescope/instrument for a given wavelength and a given set-up. We 
refer to Joos et al. (2008) for a detailed breakdown of the required 
matrix components. The Mueller matrix M(^, X) describing the ef- 
fects of the M3 mirror on the polarisation is a function of the light 
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Figure 5. This figure shows the linear polarimetry light curve with broad- 
band light curves as reference. 

The top panel shows the XRT X-ray light curve over this interval. The panel 
below that shows the r (red) and i (yellow) light curves over this time inter- 
val {i magnitudes shifted by -0.5 magnitudes). The lower two panels show 
the linear polarisation and polarisation position angle. We plot the polari- 
metric datapoints as observed, i.e. before correction for Galactic dust in- 
duced polarisation. The inset shows the circular polarimetry (four measure- 
ments of Stokes VII) performed in between linl and lin2. The horizontal 
thin red line shows the V// = level, and the blue dashed line shows the 
average of the four datapoints, V// = -0.0002, with the solid orange lines 
indicating the error 0.00075; leading to a 2cr limit 0.15%. 



incident angle 6, wavelength X and the refraction index components 
n and k of the reflecting surface of M3. For these latter values we 
use the values at the K band as derived by Witzel et al. (201 1). 

We obtained observations of the unpolarised DA white dwarf 
GJ 915 (WD 2359-434; Fossati et al. 2007) using the same instru- 
ment rotation angles as the science data and at a similar parallac- 
tic angle (60 degrees at mid of the observation vs. 50 degrees for 
the mid of the range of the science series). Our calibrator obser- 
vations were obtained just before the ISAAC science observations. 
We add to these observations from the ESO data archive of unpo- 
larised white dwarf WD 1620-391 (using identical rotator angles 
as our data), taken on 16 September 2009, with parallactic angle 
76 degrees. Using these stars and the Mueller matrices described 
above we calibrate our Stokes parameters. Using our own standard 
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Figure 6. Panels are as in Figure|5] showing the behaviour of the source dur- 
ing the first night of polarimetric observations {linl - linlO) - note the linear 
time axis. The lower two panels show the linear polarisation and polarisa- 
tion position angle, but this time we plot both the polarimetric datapoints 
as observed (red) and after correction for Galactic dust induced polarisation 
(green). Open black squares give the ISAAC Ks band polarisation degree 
and angle. Data points with zero GIP corrected polarisation have no as- 
sociated position angle measurement. We note that position angles range 
[0-180], i.e. one can add 180 degrees to green datapoints with low angle 
values (to any datapoint one can add or subtract 180 degrees) to visually 
verify they behave similarly to the trend seen in the data before GIP correc- 
tion. Dashed vertical lines indicate the centres of three bumps found in the 
optical light curve (fits are in Section |3Ti . 
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Figure 7. Panels are as in Figure |5] this time showing the behaviour of the 
source during the second night of polarimetric observations (linll - linl9). 
A 180 degree rotation of the position angle is apparent, sampled in the first 
6 datapoints. 



tion in R band of ~ 0.45% corresponds to an expected GIP contri- 
bution in the Ks band of < 0.1%, i.e. an order of magnitude smaller 
than the statistical uncertainties in the value above. We therefore 
make no correction to the Ks band Stokes parameters above. Using 
the same techniques as for the F0RS2 data we compute a bias cor- 
rected linear polarisation in Ks band of 2.0 ±0.8% with a position 
angle 10.9 + 20.9 degrees. 



star observations and those of 2008 gave consistent results within 
the errors. Because of the way we dithered our data and the low 
number of bright field stars, we can not refine our calibration using 
field stars. We then rotate the (^2, U2) values to the same orientation 
as {q, u). Because the source was fairly faint in Ks band, the errors 
on the Stokes parameters of linKl-3 are large, and we therefore 
combine them together, finding a combined q - 0.0204 + 0.0078 
and w = 0.00081 +0.0083. 

There are no field stars available for an independent GIP cal- 
ibration in Ks band, but we can estimate the likely contribution 
of GIP to the ob served values by using the approximate relation 
p{X) oc ^-18 (see lMartin & Whittetl[T99Ql) . While the value of the 
exponent has fairly large uncertainty, it is clear that a GIP polarisa- 



2.3 X-shooter spectroscopy 

After the discovery of an optical afterglow, we activated our VLT 
X-Shooter (mounted on UT2 of the VLT, Kueyen) programme 
for GRB afterglow spectroscopy (programme ID 084.A-0260, PI 
Fynbo). We obtained 4 X 600s exposures with midpoint 0.1485 days 
(1.283 X 10"^ seconds) after burst. X-Shooter is an echelle spectro- 
graph with three arms, the UVB, VIS and NIR arms, separated by 
dichroics, resulting in a wavelength coverage ~ 0.3 - lAfim. We 
used a 1 X 2 binning in the UVB and VIS arms (binning in the 
wavelength coordinate). A 5 arcsecond nod throw was used to fa- 
cilitate accurate sky subtraction, particularly important in the NIR 
arm, resulting in exposures taken in a ABBA pattern. We processed 
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Figure 8. Upper panel: P/Pmax for a range of Ry values, at the redshift of 
GRB 091018 and using host galaxy extinction Ay - \ magnitude. Lower 
panel: resulting ratio of R band and K band Hnear polarisation (in the case 
that all detected polarisation is due to dust scattering in the host galaxy). 



the data using version 1.3.0 of the ESO X-Shooter data reduction 
pipeline (Goldoni et al. 2006), using the so-called physical mode. 

Flux calibration is achieved using exposures of the flux stan- 
dard star BD+17 4708. Telluric line correction of the NIR arm 
spectrum was performed using metho ds outlined in Wiersema 
(2011 ) with the SpeXtool software ('Cushing, Va cca & Ravnej 
l2004h . using observations of the B9V star HD 16226. The resulting 
spectra are normalised. 



2.4 Optical and near-infrared photometry 

2. 4. 1 GROND afterglow photometry 

The seven channel Gamma-Ray Burst Optical and Near-Infrared 
Detector (GROND; Greiner et al. 2008) mounted on the 
ESO 2.2m telescope observed the field of GRB 091018 using 
g\r\ /', z', /, //, Ks filters beginning as soon as the source became 
visible from La Silla observatory (Filgas et al. 2009). Photometric 
calibration was achieved through observations of photometric stan- 
dard star fields. A spectral energy distribution usin g some of the 
GRON D data presented here was published earlier in lGreiner et al.l 
(1201 ih . 
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Figure 9. Emission lines of the [O II] doublet in the spectrum. The 2D spec- 
trum is shown in the top panel. Weak residuals from imperfect subtraction 
of sky emission lines can be seen as vertical stripes, also showing up as 
narrow residuals in the ID spectrum. 



2.4.2 Faulkes Telescope South afterglow photometry 

We observed the optical afterglow using the Faulkes Telescope 
South (FTS), at Siding Spring, Australia. Observations were per- 
formed at two diff'erent epochs (about 15.3 hr and 3.7 d after the 
burst event), using R and / filters. These observations are comple- 
mentary to those of GROND, filling the gap in coverage during 
daytime in Chile. All the images were cross-calibrated using the 
GROND field star measurement. 



2. 4. 3 VLT afterglow photometry 

As part of our linear and circular polarisation monitoring campaign 
described above, a large series of acquisition images were taken 
for accurate positioning of the afterglow within the aperture mask. 
Within each of these images the afterglow is detected at good signal 
to noise levels. We reduce the 14 F0RS2 and 2 ISAAC acquisition 
images using flat fields, dark- and bias fields taken the same night. 
The F0RS2 acquisition data were taken using the Rspeciai filter, the 
ISAAC data using the Ks filter. We use a sample of field stars in the 
GROND data to bring the F0RS2 acquisition photometry to the 
same homogeneous calibration. The ISAAC data were calibrated 
using four 2MASS stars in the field of view. The resulting magni- 
tudes were converted to AB magnitudes, the resulting values are 
reported in Table |3] 

In addition to the polarimetry acquisition images, we also 
photometer the acquisition image taken for the X-Shooter spec- 
troscopy, calibrating onto the GROND system. 



2. 4. 4 Gemini-South afterglow and host photometry 

We followed the late-time afterglow with GMOS-South, mounted 
on the Gemini South telescope, Chile. We observed the source at 
three epochs, in Sloan r band (CMOS filter r_G0326), see Table 
[3] We reduce the GMOS images using twilight flat fields and bias 
fields taken the same night using tasks from the Gemini packages 
in IRAF. In all three data sets the host galaxy of the burst is vis- 
ible. Photometric calibration is performed relative to the GROND 
calibration. 
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2. 4. 5 UVOT afterglow photometry 

Swiff s UltraViolet-Optical Telescope fUVOT. lRoming et alJl2005h 
began settled observations 68 seconds after the BAT trigger, begin- 
ning with a 150 s exposure in the white filter followed by a 245 s 
exposure in the u band and thereafter cycling through all 7 lentic- 
ular filters. We performed photometry and created light curves us- 
ing the tool uvotproduct provided in the Swift software. The latest 
calibration data as of July 2011 were used. We binned the data re- 
quiring a minimum detection significance of 3 cr per bin; errors are 
given at the Icr level. 



3 RESULTS 

3.1 Light Curves and spectral energy distributions 

3.1.1 X- ray light curve fits 

The Swift XRT X-ray data, as analysed through the Swift XRT cat- 
aloguqj (see Evans et al. 2009 for details), show an initial shallow 
decay phase, with power law decay index ai = 0.41^JJQg. A break 
in the light curve at 587^^^ seconds sees the decay continue with 
a steeper index 0^2 = 1.24 + 0.04. Between ~ 10"^ and 10^ sec- 
onds a broad bump is visible in the X-ray light curve. The Swift 
XRT catalogue lists the best fit solution (^^/dof = 117/133) us- 
ing two breaks in quick succession, at 2.24^q^2 x 10"^ seconds and 
2.69^Q^^ X 10^ seconds, with the first leading to a negative index 
(i.e. a brightening, but with poorly constrained index), and the last 
break leading to a final decay index of (24 = 1.59^Jj }^. It is likely 
that the behaviour between ~ 10"^ and 10^ seconds is the result of 
a system of late time flares (or some other form of short time scale 
rebrightening), in which case a more realistic alternative descrip- 
tion may be that of a single break at 4.7^^ ^ x 10"^ seconds, with 
Qf2 = 1.17 + 0.03 and (23 = 1.54^q}3. This gives a poorer ;^^ statistic 
O^^/dof = 137/135), as a cluster of datapoints interpreted as due 
to flaring/bump in this scenario are off" the best fit. The X-ray light 
curve and the fit are shown in Figure [TO] 



3.1.2 Optical light curve fits 

The optical light curves presented in this paper have their densest 
coverage through the GROND observations, as described above. 
Figures |4l[6]and|7]show the presence of bumps in the light curve, on 
top of the usual power law decay. Bumps in afterglow light curves 
are fairly common, though they require good signal to noise to be 
identified. Three distinct bumps, each lasting roughly half an hour, 
are readily visible in the first night data, and are indicated by dashed 
vertical lines in Figure |6] (see Section 13.11 these bumps are also 
detected by the PROMPT telescopes. iLaCluvze et all 2009). Simi- 
lar to the X-rays, the presence of bumps makes the light curve fits 
somewhat complex: x^ fit statistics are largely driven by datapoints 
belonging to bumps covered by GROND. To find a best possible 
way to di scern bumps and pow er law breaks, we use the method de- 
scribed in lCurran et al.l(l2007bh : data from all filters were combined 
through a free off'set fit (i.e. we assume data from all filters have 
the same temporal decay but we make no assumptions regarding 
the relative off'sets). Simulated annealing is used to fit the off'sets, 
temporal indices and break times, and a Monte Carlo analysis with 
synthetic data sets is used for error determination: the normalisation 
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error is added in quadrature to the datapoints in the resulting com- 
bined light curve. We exclude very early data (i.e. the first UVOT 
datapoint at around 100 seconds), late time UVOT data (> 10^ sec- 
onds) and the late time Gemini data, which are nearly entirely host- 
dominated. The resulting shifted datapoints are placed on an arbi- 
trarily scaled flux scale, and we fit this resulting dense light curve. 
We fit a broken power law, which gives a reduced ;^^ = 1.68 for 365 
degrees of freedom. We then add Gaussian functions one by one to 
empirically take care of the bumps in the light curve (i.e. this is not 
a physically motivated description of light curve bumps). The best 
fit is obtained using three Gaussians and a broken power law, and 
is displayed in Figure [TOl with reduced ;^^ = 1.12 (356 degrees of 
freedom). We find best fitting parameters for the broken power law 
component of ai = 0.81 + 0.01, 0^2 = 1.33 + 0.02 and a break time 
of Wak = 3.23 + 0.16 X 10"^ seconds. The three Gaussians have 
best fitting centre times 1.43 + 0.01 x 10^ 1.68 + 0.01 x 10"^ and 
2.13 + 0.06 X 10"^ seconds. Their widths as given by the standard 
deviation of the fitted Gaussians are a = 1140 + 200,370 + 80 
and 1026 + 230 seconds, respectively, though we note that the third 
bump is only partially covered by datapoints, resulting in the larger 
errors in central time and width. Adding the late time Gemini data 
to the dataset after subtracting the host brightness (the last Gemini 
epoch) does not change the fits. 

3.1.3 Spectral energy distributions and light curve interpretation 

In a study of the host galaxy dust extinction properties of sight- 
lines towards GRB afterglows using single epoch broadband spec- 
tral energy distributions (SEDs), Schady et al. (2011) report on 
a late time SED of the afterglow of GRB 091018. This SED (at 
3 X 10^ seconds, using GROND, UVOT and XRT data), is best fit 
with a broken power law, indicating the presence of a synchrotron 
cooling break in between X-ray and optical/UV wavelengths. We 
extend this analysis, fitting the X-ray - optical SEDs at three rep- 
resentative epochs: at 10^^,3 x 1 0"^ and 10^ seconds , using a fit- 
ting method described in detail in lKriihler et al.l (l201lh . We assume 
there is no change in X-ray absorption and optical extinction be- 
tween the three epochs. We fit using a broken power law (with the 
two power law indices linked with Pi - fe - 0.5) and SMC-like 
extinction law, shown by Schady et al] (l2012h to be the best fit for 
this afterglow, with the E{B - V) and A/^(H) fit simultaneously for 
all three epochs. Bands with a possible contribution from the Ly- 
man limit are excluded from the fit. The fit statistics are computed 
using C statistics in the X-rays (see Evans et al. 2009) and x^ in 
the optical. The fits are displayed in Figure [TT] The three SEDs 
are best fit ix^ = 27.04 using 31 bins, C-stat = 315.35 using 433 
bins) with the following properties (errors at 90% level confidence): 
E{B - V) = 0.024+j;;005, 7V(H) = 1.7 + 0.8 x lO^^ cm'^ and power 
law slope pi - .58 + . 07. Th ese values are in agreement with the 
results of .Schadv et al.l (120121) . We find power law break energies 
^break = 0.05:g;i3, ().()5^_^l^ and 0.025:^;^^^ keV for the SEDs at 10, 
30 and 100 ksec, respectively. 

We now compare the decay indices and spectral indices to clo- 
sure relations for some simple models. The spectral indices found 
from the SED fits imply a power law index for the electron energy 
distribution ;? = 2. 16 + 0. 14. If the blastwave propagates in a homo- 
geneous density medium, the predicted afterglow decay index for 
the regime y < y^ is 3(/7 - l)/4 = 0.87 ±0.11, and for y > y^ it 
is 1.12 + 0.11. A wind-like medium (where density decreases with 
r^) for y < Vc would imply a - (3p - 2)/4 = 1.37 + 0.11; and 
for V > Vc has the same index as for the homogeneous medium. 
The values for a homogenous medium agree well with the values 
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Figure 10. The optical light curve, in arbitrary flux units as described in 
Section [TT] is shown in blue (top), with the best fitting light curve super- 
posed: a broken power law and three Gaussian components describing the 
behaviour of flares. The red datapoints show the X-ray light curve, in units 
of count rate, fit with a double broken power law, see Section ITTI for the fit 
parameters. 
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Figure 11. X-ray - optical spectral energy distributions at three epochs (ob- 
server frame). Fits are done with a broken power law with a fixed diff'erence 
between the two power law spectral indices, and with the same X-ray ab- 
sorption and optical extinction for all three epochs. 



detennined from the light curve fits (using the two break model on 
the XRT data) for the pre-break slopes, and are consistent with the 
presence of a cooling break between X-ray and optical frequencies. 
In addition, in a wind-like medium the cooling break frequency 
should increase with time, which is not observed. 

The break time in optical and X-rays are broadly in agree- 
ment, and the SEDs (Fig[TTJ the first SED is before the break, the 
last after) show no evidence that the break is wavelength depen- 
dent (chromatic). As such, identifying this break with a jet break 
is compelling, particularly as this break is the second light curve 
break seen in the X-ray afterglow, with the first likely marking the 
end of energy injection. However, we note that the post-break de- 
cay indices (1.33 + 0.02 and 1.54 + 0.15 for optical and X-rays, re- 
spectively) are too shallow compared to predictions from the stan- 
dard fireball model, both in the case of non- sideways spreading jets 
(which predicts post-break decay indices of 1.52 + 0.11 for opti- 
cal and 2.12 + 0.11 for X-rays) and sideways spreading jets (which 
should give post-break decay indices of p). 



An alternative interpretation could be that the break is not a 
jet break, but rather a break marking the start of the "normal" after- 
glow phase, i.e. the decay indices after this break are pre-jet break 
indices, and the phase before this break is dominated by energy 
injection. In this case the post-break optical index of 1.33 + 0.02 
is consistent with the closure relations in case of a wind medium, 
which predicts 1.37 + 0.11. However, in that interpretation the op- 
tical decay index should be steeper than the X-ray one, and the 
predicted X-ray decay index of 1.12 + 0.11 is not consistent with 
the measurement. 

We note that the homogeneous analyses of large samples of 
Swift XRT afterglows have shown that a large fraction of afterglows 
with candidate jet breaks (i.e. a light curve break significantly after 
a plateau phase) produce post-break decay indices that are to o shal- 
low compared to simple model predi ctions ([Will ingale e t al.ll2007l : 
lEvans et al.l2009l . their figure 10; and lRacusin'et al.ii2009l their fig- 
ure 2). Both these papers further point out several additions to the 
models that can possibly explain some of this discrepancy, for ex- 
ample the presence of low-level continuous energy injection past 
the plateau phase, or peculiar jet structure and sightline configura- 
tions. In the following we will consider the observed break a can- 
didate jet break, and will discuss alternative options as well. 

If the break is a jet break, we can correct th e isotropic energy 
Ei^o for beaming. Using Ei^o = 3.7 X 10^^ erg taolenetskii et al. 
[2009) and making the same assumptions as in Kocevski & B utlei 
|2P08, we find a jet opening angle ^jet ~ 0.059 radians, or 3.4 de- 
grees, leading to \og(Ey) ~ 48.9. This is a fairly low value for 
Ey, but not unpre cedented: similar values are found for several 
other Swift huYsi s (iKocevski & Butlej|2008h and prc-Swift bursts 
(iFrail et aLll200lh . 



3.2 Circular polarisation 

The presence of a (weak) ordered magnetic field in addition to a 
chaotic, incoherent magnetic fie ld generated by the post - shock tur- 
bulence, has been proposed (e.g. lGranot & Koniglll20 03h to explain 
the non-detection of the swing in linear polarisation position angle 
at the time of the jet break, expected in some models (e.g. Laz- 
zati et al. 2003). In this picture, the observed polarisation may be 
dominated by the weak large scale ordered field, while emissiv- 
ity is dominated by random tangled field made by the post- shock 
turbulence. Polarisation variability occurs as a result of changes in 
the ratio of the ordered- to-random mean- squared field amplitudes 
(Granot & Koenigl 2003). This results in much weaker changes in 
the polarisation and polarisation position angle light curve around 
the jet break time. A direct test of this proposition can come from 
deep circular polarisation. Some searches have been performed at 
radio wavelengths during radio flares (Granot & Taylor 2005), but 
with fairly poor sensitivity (the best upper limit on circular polari- 
sation is 9% for GRB 991216). 

Our dataset contains four measurements of circular polarisa- 
tion in optical, R band, wavelengths, see Table [T] each with un- 
certainties of 0.15%. Under the assumption that over the interval 
that the data were obtained (~ 0.7 hours, see Table [B there is no 
change in circular polarisation, we can combine these together. We 
find a combined value of the circular polarisation Stokes parameter 
of V = V/I = -0.00020 + 0.00075, which leads to a formal Icr up- 
per limit of Pcirc < 0.15%: the deepest limit on circular polarisation 
of a GRB afterglow to date. Figures |4] and [6] show that during the 
circular polarimetry epoch, the optical light curve shows a low am- 
plitude bump. As such, the limit on the circular polarisation can be 
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seen as a limit on the circular polarisation of the light of the bump 
plus that of the underlying afterglow. 



3.3 Interstellar polarisation within the host galaxy 

The variability of the detected linear polarisation indicates that the 
majority of detected polarisation can be associated with the after- 
glow. However, scattering of afterglow light onto dust grains within 
the host galaxy (HGIP (host galaxy interstellar polarisation), fol- 
lowing terminology of Gorosabel et al. 2010) can induce notice- 
able linear polarisation, depending on the dust scattering geometry 
and grain size distribution. This can affect a ttempts to interpret po- 
larisation behaviour around jet breaks (e.g. iLazzati et al.ll2003h or 
models using the absolute level of polarisation (e.g. Gruzinov & 
Waxman 1999; Gruzinov 1999). For example, the zero polarisation 
seen in e.g. linl could in principle be an unfortunate effect of HGIP 
on an intrinsically non-zero polarised afterglow. 

We can utilize the reasonable assumption that over the opti- 
cal range, the intrinsic afterglow linear polarisation will be wave- 
length independent as it is synchotron emission, whereas HGIP 
will be strongly wavelength dependent. In the Galactic ISM, the 
polarisation broadly follows an (empirical) relation known as the 
Serkowski law: 



^lin = ^lin,max CXp UCXn^ \~^]) ' 



In this relation, Piin,max is the maximum induced linear polarisation 
at wavelength Aj^^x- The wavelength /Imax traces the size (distribu- 
tion) of the dust grains responsible for the observed polarisation, 
and as such is closely linked to Ry = Av/E(B - V). Klose et al. 
(12004) considered the effect of redshift in the Serkowski law, which 
results in a substitution /l^ax -^ (1 + z)A^^l in the equ ation above, 
and the relation A^^l = Rv/5.5. iKlose et alJ (l2004l) derive the 
observed linear polarisation P\m/P\m,max in standard photometric 
broadband filters of an intrinsically unpolarised afterglow due to 
the presence of dus t in the host , as a f unction of redshift z and Ry. 
Using the results in lKlose et al.l (120041) . we generate predictions for 
observed Piin/^iin,max in several broadband filters as a function of 
Ry, using z = 0.971, shown in Figure [8] From these, a predicted 
ratio of Piin,R/^iin,K is produced as a function of Ry, plotted in the 
lower panel of Figure [8] Multi- wavelength linear pola rimetry has 
been possible for only a small number bursts: 020813 (JBarth et al.l 

t3; Lazzati et al. 2004), 021004 (Laz zati et al.l2003h and 030329 
)se et alJl2004l : lGreiner et al.l l2003). In these cases, no signifi- 
cant evidence for HGIP has been observed. 

Our Ks band polarimetry is approximately simultaneous to 
the R band polarimetry epoch lin9 (to be exact, lin9 is simul- 
taneous with linK3). The (GIP corrected) Stokes parameters for 
lin9 and linK are very similar, though we caution that the data 
were taken during what appears to be a bump in the polarisa- 
tion light curve (Fig. [6]). The resulting position angles are iden- 
tical within errors, and we find Piin,R/^iin,K = 0.72 + 0.45. This, 
combined with the strong detected variation of both polarisation 
degree and angle, demonstrates that (within errors) there is no in- 
dicat ion of substantial HGIP contribution to the observed polarisa- 
tion. iMartin & Angell (Il976h demonstrate that the ISM induces a 
small degree of circular polarisation even in unpolarised sources. 
However, the levels of induced circular polarisation are an order of 
magnitude or more below the limit we set in Section 13.21 consid- 
ering the measured Galactic GIP and the low extinction in the host 
galaxy (Section [3T|). 

In the following we therefore proceed from the assumption 



that the GIP-corrected polarisation is all intrinsic to the source. In a 
future paper (Paper 2; Wiersema et al. 2012 in prep.) we will further 
quantify the implication of this measurement (and those of other 
GRB afterglows) on the dust and gas properties in GRB sightlines. 
A possibility that we will also address in Paper 2 is the effect of 
dust destruction by the GRB and afterglow high energy radiation, 
which may make the effect of HGIP time-dependent. 



3.4 Linear polarisation light curve 

Our R band linear polarimetric monitoring took place over three 
observing nights, as described in Section IT. 1.1 1 The behaviour over 
the first two nights is illustrated in Figures |6]and|7] shown in linear 
time coordinate. On the third night only one deep, final, measure- 
ment was taken. 

In the first night (see Figure |6j, the source starts off with very 
low (consistent with zero, within the errors dominated by the un- 
certainty in the GIP Stokes parameters) polarisation, the lowest ob- 
served for a GRB afterglow to date. Following that, a broad, low 
amplitude bump may be present in the polarimetric curve. GRBs 
030329 and 021004 showed a correlation between linear polar isa- 
tion behaviour and light curve bumps (^e.g. iGreiner et al.ll2003h . In 
Figure [6] we indicate the times of the three bumps we find in our 
data. Some variability may be present related to these bumps, but 
the statistical significance is very low. After the first eight data- 
points obtained at the start of the night, observations were resumed 
at the end of the night, showing a decline from a higher polarisation 
level (Table 2), confirmed by the ISAAC polarimetry taken quasi- 
simultaneous to lin9. Over the first night, some low level variation 
in position angle may be present, but the uncertainties are large be- 
cause of the low linear polarisation: overall, the first night data are 
consistent with a constant polarisation angle. 

The second night data (see Figure IT} displays short time scale 
variability: a polarisation bump, peaking at > 3% polarisation, ac- 
companied by a clear rotation of the position angle. Following this, 
the end of the night shows a slow rise in polarisation, accompa- 
nied by further position angle variation. By the third night, the af- 
terglow had become too faint for detailed sampling, and one deep 
datapoint, linlO, was obtained (Table 2). This last datapoint clearly 
shows low-level contribution from the host galaxy in the received 
light. We will attempt to interpret the polarimetric data in Section 

m 



3.5 Environment of the burst 

In a forthcoming paper on this dataset (Paper 2) we will analyse 
the GRB environment in detail. In this section we will derive some 
properties most relevant to the afterglow physics, and to show that 
this burst displays no indications of being unusual. 

In the late time images taken with Gemini (see Table [3]) the 
host galaxy can clearly be seen as an extended structure, see Figure 
[T2I visibly extended in East- West direction. An isophote fit using 
elliptical isophotes (using the ellipse package within IRAF) gave 
a best fit for an ellipse orientation of 89 + 4 degrees. We deter- 
mined the position of the GRB within its host by performing image 
subtraction of the three Gemini epochs using a modified version 
of the ISIS2 code (Alard & Lupton 1998). We detected a clear af- 
terglow residual in a subtraction of the first and third epoch, and 
between the first and second epochs. No residual was detected in 
a subtraction between the second and third epochs. We found the 
GRB location to be offset by approximately 5.8 kiloparsec from 
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the centre of the galaxy (as found from the eUiptical isophote fits), 
see Fig [121 This offset, as well as the observed host magnitude, 
r = 23.4 (AB; see Table[3]), is well within the observed distribution 
for long GRBs at this redshift fe.^.lBloom. Kulkami & Djorgovski 
l2002l : ISavadio. Glazebrook & LeBorgnell2009h . 

The X-shooter afterglow spectra show very strong absorption 
lines (Fig (Ts]) composed of at least three discrete velocity com- 
ponents, spanning ~ 300 km s~^ No lines from intervening sys- 
tems are detected. At the GRB redshift we detect both common 
resonance lines and a small number of transitions arising from ex- 
cited finestructure and metastable levels (from Fell and possibly 
Nill), though the lines from these excited transitions are very weak 
in our spectrum. Time resolved spectroscopy of GRB afterglows 
has shown that these excited transitions are excited through in- 
direct UV-pumping (fl uorescence) by the afterglow photons (e.g. 
IVreeswiik et al.ll2007h . This unambiguously identifies the redshift 
z = 0.971 to be that of the GRB. 

The X-shooter spectra were taken with the slit position an- 
gle aligned with the parallactic angle, set at -90.8° at the start 
of observations, and therefore oriented broadly along the host 
galaxy. In addition to absorption lines, we therefore also detect 
nebular emission lines in the spectrum (e.g. Ha, [Oil], [OIII]), 
though their signal to noise is low as the afterglow strongly dom- 
inates the light (see Table [3]). The detection of the resolved [Oil] 
.1.13726, 3729 doublet is shown in Figure|9] Its flux of ~ 6 x lO'^^ 
erg s~^ cm~^ can be used to derive an estimate of the starforma- 
tion rate: using the equation from Kennicutt (1998) we find a 
SFR([OII]) ~ 4 Solar masses per year. If we instead use the 
empirical conv ersion from [O II] line luminosity to star formation 
rate derived bv lSavadio. Glazebrook & LeBorgnel(l2009h for GRB 
host galaxies, we find SFR([OII]) ~ 1.6 Solar masses per year 
(note that no extinction correction was applied to the [Oil] line 
luminosity). This SFR([OII]) is fairly typical for long GRB hosts 
(I Savaglio. Glazebrook & LeBorgnell200 9). In conclusion, the after- 
glow spectrum, host magnitude and morphology, and star formation 
rate are all in line with expectations for a normal long GRB host 
galaxy. 

In Paper 2 we will further analyse the spectrum, and exploit 
the rare opportunity where we can probe (interstellar) gas and dust 
in the host galaxy using five methods simultaneously: through after- 
glow absorption lines (including excited finestructure lines), nebu- 
lar emission lines, afterglow SEDs, wavelength dependent linear 
polarimetry and host imaging. 



4 DISCUSSION 

4.1 GRB 091018 compared to GRBs 030329, 021004, 020813, 
020405 and 090102 

It is instructive to compare the properties of the polarisation of 
the afterglow of GRB 091018 with those of the best studied pre- 
vious cases: the pve-Swift bursts 030329 (z = 0.1687), 021004 
(z = 2.323), 020813 (z = 1.254) and 020405 (z = 0.690). Each has 
multiple measurements over more than one night. We also compare 
to Swift burst 090102 (z = 1.547), which has only one single mea- 
surement, which probes the poorly explored early-time afterglow. 

Both 021004 and 030329 h ave well s tudied bum ps in their 
afterglow lig ht curv es (021004: ' Bersier e t al. 2003bl; IFox et al.l 
1200 3: Hollan d etal.|[200 3: Mirabal et al. 2003; Pandev et al. 2003. 
jm et al.l2004,:.Burenin et al...2003.:.L ipkin et al. .2004, : 
an i2003l : 




030329jBIoom 



iMatheson et"ani2003l : [Sato et al.|[2003l : lUemura et al. 2003]), fre- 
quently attributed to late-time energy injection (late time activity 



Figure 12. The host galaxy of GRB 091018, from the second epoch Gemini 
GMOS r band data. The dashed tickmarks mark the position of the GRB as 
determined from PSF-matched image subtraction of Gemini epochs 1 and 

2. 



of the central engine). Both sources are pre- Swift bursts and there- 
fore have relatively poorly sampled X-ray light curves compared 
to bursts in the Swift era, but well sampled optical and radio after- 
glows. GRB 030329 has the best sampled polarimetric light curve 
of any burst to date, with polarisation measurements extending to a 
month after burst (.Greiner et al...2003.) . The data show some corre- 
lation between the degree of linear polarisation Pun and changes in 
the decay index of the afterglow (bumps), though Pun and the light 
curve do not change strictly simultaneously (Greiner et al. 2003|). 
Similarly, it appears that not every bump is associated with a polar- 
isation change. Variability in both p osition angle and polarisation 
degree is detected for 021004 (e.g. IRoI et al.ll2003l: iLazzati et aJ] 
2003), but an association with flare behaviour is not easy to make, 
in part because the polarimetric coverage is less dense than for 
030329. In contrast to 021004 and 030329, the afterglow of 020813 
is smooth, and therefore this afterglow may be better suited to study 
global jet properties. Similarly to 021004 and 030329, variability 
in polarisation is detected in 020813. GRB 091018 is the first Swift 
burst to join the small group of afterglows with observed polarimet- 
ric variability. Bumps are observed in the broadband optical light 
curves of this source, but their amplitude is considerably less than 
those seen in the afterglows of GRBs 021004 and 030329. If the 
bump amplitude correlates with polarisation amplitude, we may 
expect the bumps to have a smaller eff'ect on the observed over- 
all polarisation, and be in a better position to probe the properties 
of the jet as a whole. 

GRB 020405 is of interest mainly because of the short time 
scale variability with high polarisation amplitu de (nearly 10% po- 
larisation) that thi s source appears to show (iCovino et al.l l2003l : 
iBersier et al.l 2003), though ne ar simultaneous measurements de- 
scribed in Masetti et al.l (l2003h do not confirm this behaviour. No 
jet break was seen in the afterglow light curve, and the polarisation 
curve could not be reconciled with standard models. 

Jetted fireball models predict a clear change in polarisation 
properties around the jet break time. The behaviour of Pun and 
the polarisation position angle depends sensitively on the angle be- 
tween the jet axis and the observer, the internal structure of the jet, 
the magnetic field properties within the jet and whether or not side- 
ways spread ing of the jet occurs at around the jet break time (e.g. 
iRossi et alj 2004 and references therein). Attempts to fit these mod- 
els to afterglow polarisation data have been a little disappointing: 
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Figure 13. The UVB arm spectrum from X-shooter, smoothed with a 3 pixel boxcar for presentation purposes. Some commonly observed (resonance) 
absorption lines are indicated with red, solid, vertical lines and are labelled. The dash-dotted, blue, vertical lines mark the positions of excited finestructure 
transitions of Fe II, which are occasionally observed as highly variable lines in GRB afterglow spectra. The dotted, green, vertical lines mark the positions of 
excited finestructure transitions of Ni II. 



in 030329 and 02100 4 the bumps in the l ight curve may have influ- 
enced the fits (see e.g. Lazzat i et al.l2003h . whereas the polarimetric 
sampling of 020813 was rather too low to favour a single jet model, 
though several models could be excluded (ILazzati et al.ll2004h . The 
afterglow of 020405 unfortunately shows no jet break within the 
interval of polarimetric monitoring. The very well sampled 030329 
has the relative disadvantage of a bright supernova contributing at 
late times, in addition to the bumps at early times, and determin- 
ing a unique jet break time proved com plicated (with best model 
likely requiring two jet breaks, see e.g. IVan der Horst et al.1l2008l 
and references therein). Our data of 091018 may therefore be one 
of the first cases where a jet break can be probed in some detail 
through polarimetry, if the light curve break discussed in Section 
13. H is indeed a jet break. 

Our polarimetric followup starts at restframe time 0.067 days 
(5782 seconds) after burst trigger, considerably earlier than the ear- 
liest datapoints of 030329, 021004 and 020813. It is clear that po- 
larimetry very early after burst requires automated/robotic proce- 
dures. Using such an approach, a very ear ly afterglow polari metric 
detection of GRB 090102 was obtained bv lSteele et all ([2009 \ who 
secure a detection of Pun = 10.2 + 1.3% at restframe time 75 sec- 
onds after burst, using the robotic Liverpool Telescope. Such a high 
degree of polarisation, measured when the reverse shock compo- 
nent was still bright, is consistent with the presence of large-scale 
ordered magnetic fields. Although our data on 091018 are taken 
considerably later in time than those on 090102, the superb sensi- 
tivity of the VLT may allow us to continue to probe the magnetic 
field structure at late time through our high-precision linear and 
circular polarimetry. 



4.2 GRB 091018 and models 

4. 2. 1 Reverse shock 

The start of our polarimetric monitoring was sufficiently late that 
a contribution of reverse shock emission to the received light is 
likely negligible. While this precludes study of some magnetisa- 



tion properties fe.g. lGranot & Tavlorll2005h . this makes a compari- 
son of the data with models easier, as only a forward shock model 
needs to be considered. The low linear polarisation at the start 
of our monitoring appears to confirm this: emission from reverse 
shocks is expected to be considerably polarised, as the shock trav- 
els back into magnetize d shocked material (e.g. iGranot & Tavlorl 
l2005l : ISteeleetal.ll2009h . 



4.2.2 Jet break 

We compa r e our data with the model series computed by 
iRossi et al.l (120041) . hereafter referred to as R04, who derive the 
expected linear polarisation curves (including position angle) for 
a grid of different jet model parameters, assuming there is no ad- 
ditional coherent component of the magnetic field present, and as- 
suming HGIP is negligible. In particular, R04 consider the cases 
of homogeneous jets, structured jets and Gaussian jets, each with 
a grid of different viewing angles (in the case of homogenous jets 
this is the angle the line of sight makes with the jet axis), and for a 
range of wavelengths and physical parameters. They further distin- 
guish jet breaks with and without sideways spreading. 

We can briefly summarize the main diff'erences between the 
various models as found by R04, and refer for details and explana- 
tions to that paper (see also Figure 18 of R04): 
(i): Homogeneous jets show two peaks in the polarisation curves, 
with the second peak always reaching a higher polarisation than the 
first. Structured and Gaussian jets show one peak only, with peak 
time for a Gaussian jet significantly later than for a structured jet, 
which has its peak around the time of the break. The peaks of the 
structured and Gaussian jets are much wider than the ones for ho- 
mogeneous jets (i.e. rise time and decay time scales are long); 
(Uy. The position angle for a homogeneous jet changes by 90° in 
between the two peaks in the polarisation curves: the two peaks 
each have constant position angle, but are 90° diff'erent. Structured 
and Gaussian jets have constant position angles throughout; 
(iUy. Structured jets show non-negligible polarisation at early times. 
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Figure 14. All GIP-corrected linear polarisation data plotted as a function 
of the time since burst divided by the time of the break in the optical light 
curve (i.e. Veak = 3.23 x 10^ seconds). Overplotted in blue is the model 
(from lRossi et al.l2004r) for a homogeneous jet with a jet break at tbreak that 
does not experience sideways expansion, with viewing angle ^obs = 0-2^jet- 
The model has been shifted sHghtly (by dividing the model f/Zjet by 1.15). 
The lower panel shows the polarisation position angle. The dotted line for 
^/^jet < 1.5 shows the average position angle (6 degrees). The dotted line 
for t/tjQt > 1.5 is set at 96 degrees: the average value before this time plus 
90 degrees, to show the predicted 90 degree position angle swing. The inset 
shows the rapid polarisation rotation at /^/Veak -3.16 (vertical dotted line 
in inset), with the horizontal dotted line at a position angle of 96 degrees. 



whereas homogenous jets and Gaussian jets have zero polarisation 
at early times; 

(iv): The shape of the bumps in the different jet structures is differ- 
ent; 

(v): Sideways expanding jets show lower polarisation values than 
non sideways expanding jets (all other parameters being equal), 
with wider bump(s) in the polarisation curve. 

We start by evaluating the data in time coordinates normalised 
by the time of the break seen in the optical data (see Section [3Jl . 
It is clear that some polarisation behaviour is observed occurring 
close to ^/^break ~ 1, as cxpcctcd for jet breaks. In the case that 
this break is a jet break, we can then directly compare the curves 
with the ones in R04. Fig. [14] shows a plot in these units. Com- 
paring to the points / - v listed above, there appears to be evidence 
that: a model with constant position angle throughout will not pro- 
vide a good fit; early time polarisation is consistent with zero (or 
at least is lower than predictions for structured jets); a single bump 
model with a shape as predicted by the models in R04 seems in- 
consistent with the data. On Fig. [141 we therefore plot a model for 
a homogeneous, non-sideways expanding jet with viewing angle 
^obs = 0.2^jet (from Figure 8 in R04). We shift this curve slightly 
in time coordinate by dividing t/tj^t by 1.15 to get a better visual 
match. The first bump in the polarisation curve appears to be well 
reproduced in both shape and amplitude, including the zero polar- 
isation at early times. This model predicts a constant polarisation 
position angle over the first bump. We indicate in Fig. [14] (dotted 
line) the average of the position angle at ?/Wak ^ 1-4. As required 
by the R04 model, the position angle is constant over this bump. 

The second bump predicted by the homogeneous jet model 
should have a larger amplitude than the first, a position angle 90° 
offset from the first bump, and a morphology opposite to the first (in 
coordinates used in Fig[T4l). Some of the data points indeed have a 
higher polarisation than seen in the first bump, but the morphology 
of the bump as predicted by the model is not observed. We examine 
the behaviour of the position angle by adding 90° to the position an- 
gle for ^/^break ^ 1-4. In Fig.[T4]this is the second horizontal dotted 
line. It is clearly seen that the value of the position angle is overall 
higher than in the first bump, but is not constant as the model pre- 



dicts, showing considerable variability on a short time scale (possi- 
bly varying around the predicted position angle), including a rota- 
tion of the position angle visible in the interval ^/ Wak ~ 3.0-3.3. 
The polarisation curve shows, compared to the homogeneous jet 
model, a deep short dip, in the interval ^/^ak ~ 3 - 4. This dip 
shows a very steep decay in the first few datapoints, simultaneous 
with the observed rotation in position angle. This behaviour is not 
predicted by any jet break model. 

We could consider the possibility that there is an additional 
polarisation curve component (hereafter simply referred to as the 
extra component) on top of a slower behaviour as predicted by jet 
break models. A third smooth, fairly narrow, polarised bump peak- 
ing at around ^/^break ~ 3 - 4 could be present, having a position 
angle significantly offset from the 6 = 96° of the second bump 
from the homogenous jet model. The sum of this component and 
that of the (second bump of the homogenous) jet model can give 
rise to most of the observations: the steep dip, where the contri- 
butions of the two components nearly cancel due to differing po- 
sition angles, and the rotation in position angle, where the com- 
ponents are at similar amplitude but very different position angle. 
However, this is a rather degenerate problem, as the amplitude as 
a function of time (the shape) of the extra component and the po- 
sition angle as a function of time are all unknown. A component 
with a fixed position angle but variable amplitude (i.e. a compo- 
nent much like the polarisation peaks from a jet break) can ap- 
proximate the data (with a position angle ~ 70 degrees offset from 
the constant 96°, and peaking around ^/^break ~ 3.5), but a com- 
ponent with nearly constant amplitude but position angle varying 
on very short t ime s cales can provide a reasonable approximation 
too. IWu et alj (120051) evaluate the expected polarisation curves of 
afterglows from tw o-component jets (as proposed for e.g. 030329, 
iBerger et al1l2003h . which show, as expected, a large variety of be- 
haviour, depending mostly on viewing angle and contrast between 
the narrow and wide jet components. However, the observed vari- 
ability in the data of 091018 around ^/^break ~ 3 - 4 is much more 
rapid than predicted by the models of IWu et al.l (1200 5). At the cur- 
rent data quality, and including the uncertainty in and applicability 
of the underlying model (homogeneous jet model from R04), we 
are not in a position to quantify the behaviour of this putative com- 
ponent further. As the shape of the bump and time variability of the 
position angle are unclear, it is difficult to speculate on the possible 
origin of such an additional component. The feature occurs after 
the expected time of the peak of the second bump in the R04 model 
and no such feature is visible at ? < ^break- As such, it may be related 
to a discrete bright region (patch) on the jet that comes into view 
just after the second bump peaks. In this scenario we may expect a 
constant position angle of the component produced by the patch (if 
the field in the patch is coherent), with a bump shape reflecting the 
brightness profile of the bright patch. However, other explanations, 
such as microlensing, two-component jets or others, are difficult to 
rule out with the polarisation data alone. 

The time span possibly most affected by the extra component, 
t/hreak ~ 3 - 4, corrcsponds to 9.7 - 12.9 x 10"^ seconds. In the op- 
tical light curves no significant slow variability at around that time 
is obvious. The XRT data show a possible broad, slow softening 
around this time, but our SED fits do not point to strong spectral 
variation at this time. The shallow post-break decay indices and/or 
the flares observed in the optical light curve (Section 13. U may be 
related to the observed polarisation behaviour, i.e. be linked to the 
observed extra component, but further modelling will be required 
to test this. 

There is at least one other source that shows a variability in the 
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polarisation curve on short time scales (intranight) and unrelated to 
a jet break, GRB 020405. However, in that source only a single 
datapoint showed the deviation from a constant polarisation degree 
(IBersier et al. 2003), and the position angle was not seen to vary 
over the available polarisation measurements. No features resem- 
bling jet break bumps were detected in that source, consistent with 
the absence of a jet break in the light curves. The possible detection 
of a variation at fairly short time scales, with possibly a constant po- 
sition angle and likely unrelated to a jet break in GRB 091018 (the 
extra component in the polarisation curve) may be an indication 
that this type of behaviour is indeed real. 

The conclusion from the above discussion is that none of the 
simple models of polarisation behaviour of GRBs directly repro- 
duces the observed behaviour. Of the available models, the homo- 
geneous jet break model comes closest, but requires at least one 
additional component, the nature of which is currently unknown. 



4. 2. 3 Ordered fields 

The circular polarisation upper limit derived in Section 13.21 
can be used to probe large scale ordered fields, which have 
been proposed by several authors to explain various phenom- 
ena, such a s the high linear polarisation at very early times 
JSteele et al.l l2009h . or the lack of obvious polarisation posi- 
tion angle shifts ar ound the jet break ti me for e.g. GRBs 
030329 and 021004 (iGranot & Koniglll2003h . Models of the ex- 
pected circular polarisation in optical wavelengths for an ex- 
ternal shock (Matsumiva &Ioka 200 3; Sagiv, Wax man & Loebl 
|2004 iToma. loka & Nakamural l2008h generally find that after- 
glows should show very low circular polarisation in the ab- 
sence of strong large scale ordered fields. Our detection limit of 
^circ < .15% is in full agreement with these predictions (see e.g. 



iToma. loka & Nakamura 200 8), and excludes strong ordered field 
components. We note the timing of our circular polarimetry: our 
data were taken when the reverse shock emission is likely negli- 
gible, and when the linear polarisation was still very weak. It is 
clear that it is not easy to significantly improve on the depth of this 
limit in optical wavelengths, but it would be interesting to search 
for circular polarisation in other bursts, to establish whether this is 
a common feature, and to search for circular polarisation at early 
times when reverse shock emission is still present. 



4.2.4 Energy injection 

An alternative explanation of the observed break in the light curves 
is that this break marks the end of (continuous) energy injection, 
though the post-break decay indices are not as expected from the 
closure relations (Section lTll . The occurrence of significant polari- 
metric variability at around f/Wak ~ 1 indicates that if this scenario 
is correct, the end of energy injection must be accompanied by a ge- 
ometric change, as the polarisation position angle stays constant. In 
a study of GRB 021004, Biomsson, Gudmundsson & Johannesson 
(I2OO4I) model the strong rebrightenings visible in the afterglow 
light curve with discrete, powerful, energy injections, and demon- 
strate that the light curves are fairly well reproduced this way. 
They demonstrate that their model shows some correlation with 
the polarimetr ic data, i.e. that polar isation is related to bumps. For 
GRB 030329, iGreiner et al.l (l2003h come to a similar conclusion 
from their data. We see at least three low amplitude bumps in the 
optical light curve, at times given in Section [TT] but a correlation 
with polarimetric behaviour is not clearly visible, largely because 



of the larger uncertainties in our polarisation values compared to 
the much brighter afterglow of 030329. Similarly, a possible re- 
lation between the extra component in the polarisation curve and 
the flaring around the jet break time is unclear. Future (n umerical) 
models of light cu rve bumps/flares (e.g . lviasis et al.l201lh and high 
quality data (e.g. iKriihler et al.ll2009h may provide the possibili- 
ties to compute expected polarisation signals from the simulations. 
These may then be used to understand the connection between the 
physical parameters behind these bumps and the behaviour in po- 
larised light. 



5 CONCLUSIONS 

We presented a large polarimetry dataset of the afterglow of 
GRB 091018, and add to that an extensive multicolour photometric 
dataset and medium resolution spectroscopy. We find an abundance 
of absorption lines in the spectrum of the afterglow as well as a few 
emission lines, from which we find a redshift z - 0.971. The X- 
ray light curve shows at least two breaks, as well as a complex 
of flaring activity. The densely sampled optical light curve shows 
several flares as well, and through a fit using a broken power law 
and three Gaussian components (to describe the flares), we find an 
achromatic break at fbreak ~ 3.2 x 10"^ seconds. Our polarisation 
data cover this break, allowing us to probe its nature using linear 
polarimetry. Our polarimetry results can be summarised as: 

• In contrast to GRBs 030329 and 021004, there are no strong 
rebrightenings that confuse our view of the polarisation behaviour. 

• Our deep circular polarimetry at f ~ 0.15 days after burst 
shows no circular polarisation, with Icr limit Pcirc < 0.15%. 

• The quasi- simultaneous R and Ks band polarimetry shows no 
strong indications of a considerable contribution of polarisation 
from the ISM in the host galaxy. 

• The linear polarimetry shows a constant position angle for 
tlhxQsk < 1-4, and a considerably higher, but strongly variable, value 
thereafter. 

• At least two peaks are present in the polarisation curve. The 
first of these, as well as the constant position angle, is well repro- 
duced by a model of a jet break using a homogenous, non-sideways 
spreading jet with small viewing angle (from R04). 

• If the observed break is a jet break, and if the model of a 
homogenous, non-sideways spreading jet is correct, an additional 
component is required to reproduce the observed polarisation curve 
and position angles. This may be a short duration bump with con- 
stant position angle, off'set from the position angle of the underlying 
jet break component. 

From the above it is clear that the models for polarisation be- 
haviour around the time of jet break (as in R04) can not, in their 
current form and on their own, explain the full behaviour of the po- 
larisation curve of 091018. Further modelling eff'orts, in particular 
into the eff'ects of rebrightenings / bumps and continuous energy 
injection on the observed polarisation, are vital to exploit the full 
potential of polarimetry as probe of afterglow physics. 

These results are the first of their kind for Swift GRBs, and in- 
deed GRB 091018 can be considered a "normal" GRB, with none 
of its prompt emission or afterglow properties particularly remark- 
able, in contrast to some of the polarimetrically studied pre-5w//f 
bursts like 030329. This is in part because our selection criteria of 
sources to study with polarimetry were only that: the source had to 
be a Swift burst with a detected afterglow in both XRT and UVOT 
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data. This avoids a bias towards the brightest or slowest decay- 
ing optical afterglows, and indeed it is clear from Section |3J1 that 
GRB 091018 is one of the faintest afterglows with multi-night po- 
larimetry. Further, similar, observations of (Swift) bursts with dif- 
ferent micro- and macro-physical parameters (in particular ener- 
getics, jet opening angle, density n and density gradient k) will be 
required to establish whether or not GRB 091018 is representative 
of the sample of GRB afterglows as a whole. 

Unfortunately, no long wavelength (radio, (sub-)mm) obser- 
vations were available for GRB 091018. A deep polarimetric cam- 
paign on an afterglow, with radio and (sub-)mm light curve data in 
addition to optical and X-ray light curves would be particularly in- 
teresting, since with the resulting broadband SEDs (with sensitivity 
redwards of the synchrotron peak frequency y^) we would be able 
to fit for parameters that are particularly relevant for polarimetric 
interpretations, e.g. €e and 65. 
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Table 3. Log of our photometry. All magnitudes are AB magnitudes. These 
magnitudes are not corrected for the Galactic foreground extinction. ^^^: cal- 
ibrated onto Sloan r using GROND. 
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0.15551 


115 


19.15 + 0.02 




0.15762 


115 


19.14 + 0.02 




0.15985 


115 


19.15 + 0.02 




0.16216 


115 


19.16 + 0.02 




0.16429 


115 


19.17 + 0.02 




0.16648 


115 


19.17 + 0.02 




0.16867 


115 


19.16 + 0.02 




0.17096 


115 


19.17 + 0.02 




0.17318 


115 


19.22 + 0.02 




0.17541 


115 


19.23 + 0.02 




0.17765 


115 


19.26 + 0.02 




0.17998 


115 


19.32 + 0.02 




0.18214 


115 


19.33+0.02 




0.18437 


115 


19.34 + 0.02 




0.18660 


115 


19.34 + 0.02 




0.18891 


115 


19.35 + 0.02 




0.19115 


115 


19.34 + 0.02 




0.19333 


115 


19.36 + 0.02 




0.19552 


115 


19.38 + 0.02 




0.19781 


115 


19.39 + 0.02 




0.19996 


115 


19.42 + 0.02 




0.20214 


115 


19.46 + 0.02 




0.20433 


115 


19.47 + 0.02 




0.20667 


115 


19.48 + 0.02 




0.20886 


115 


19.48 + 0.02 




0.21105 


115 


19.51+0.02 




0.21324 


115 


19.54 + 0.02 




0.21556 


115 


19.54 + 0.02 




0.21767 


115 


19.53+0.02 




0.21985 


115 


19.54 + 0.02 




0.22204 


115 


19.54 + 0.02 




0.22433 


115 


19.53 + 0.02 




0.22648 


115 


19.51+0.02 




0.22867 


115 


19.52 + 0.02 




0.23086 


115 


19.54 + 0.02 




0.23313 


115 


19.54 + 0.02 




0.23526 


115 


19.54 + 0.02 




0.23749 


115 


19.54 + 0.02 




0.23977 


115 


19.54 + 0.02 




0.34860 


675 


20.01 + 0.02 




0.44872 


688 


20.21 + 0.02 




1.20917 


675 


21.76 + 0.04 




1.30936 


688 


21.82 + 0.04 




1.49485 


688 


22.02 + 0.05 




2.22673 


1727 


22.53 + 0.04 




3.15788 


1722 


22.94 + 0.07 


GROND, r 


0.12693 


115 


18.87 + 0.03 




0.12906 


115 


18.86 + 0.02 




0.13125 


115 


18.85 + 0.02 




0.13344 


115 


18.87 + 0.02 




0.13575 


115 


18.88 + 0.02 




0.13793 


115 


18.91+0.02 




0.14012 


115 


18.91+0.02 




0.14231 


115 


18.92 + 0.02 




0.14457 


115 


18.93 + 0.02 




0.14670 


115 


18.93+0.02 




0.14889 


115 


18.94 + 0.02 




0.15109 


115 


18.94 + 0.02 




0.15330 


115 


18.94 + 0.02 




0.15551 


115 


18.96 + 0.02 




0.15762 


115 


18.97 + 0.02 




0.15985 


115 


18.98 + 0.02 




0.16216 


115 


18.98 + 0.02 
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0.16429 


115 


18.99 + 0.02 




0.16648 


115 


18.98 + 0.02 




0.16867 


115 


18.97 + 0.02 




0.17096 


115 


18.99 + 0.02 




0.17318 


115 


19.03 + 0.02 




0.17541 


115 


19.06 + 0.02 




0.17765 


115 


19.08 + 0.02 




0.17998 


115 


19.11 +0.02 




0.18214 


115 


19.14 + 0.02 




0.18437 


115 


19.15 + 0.02 




0.18660 


115 


19.17 + 0.02 




0.18891 


115 


19.15 + 0.02 




0.19115 


115 


19.17 + 0.02 




0.19333 


115 


19.17 + 0.02 




0.19552 


115 


19.18 + 0.02 




0.19781 


115 


19.21 + 0.02 




0.19996 


115 


19.24 + 0.02 




0.20214 


115 


19.25 + 0.02 




0.20433 


115 


19.27 + 0.02 




0.20667 


115 


19.30 + 0.02 




0.20886 


115 


19.32 + 0.02 




0.21105 


115 


19.33 + 0.02 




0.21324 


115 


19.33 + 0.02 




0.21556 


115 


19.37 + 0.02 




0.21767 


115 


19.36 + 0.02 




0.21985 


115 


19.37 + 0.02 




0.22204 


115 


19.36 + 0.02 




0.22433 


115 


19.33+0.02 




0.22648 


115 


19.34 + 0.02 




0.22867 


115 


19.35 + 0.02 




0.23086 


115 


19.35 + 0.02 




0.23313 


115 


19.34 + 0.02 




0.23526 


115 


19.35 + 0.02 




0.23749 


115 


19.35 + 0.02 




0.23977 


115 


19.37 + 0.02 




0.34860 


675 


19.82 + 0.02 




0.44872 


688 


20.01 + 0.02 




1.20917 


675 


21.59 + 0.03 




1.30936 


688 


21.68 + 0.03 




1.49485 


688 


21.88 + 0.05 




2.22673 


1727 


22.34 + 0.04 




3.15788 


1722 


22.81+0.10 


GROND, i 


0.12693 


115 


18.71+0.04 




0.12906 


115 


18.71+0.03 




0.13125 


115 


18.72 + 0.02 




0.13344 


115 


18.75 + 0.03 




0.13575 


115 


18.76 + 0.03 




0.13793 


115 


18.79 + 0.02 




0.14012 


115 


18.79 + 0.02 




0.14231 


115 


18.78 + 0.02 




0.14457 


115 


18.81+0.03 




0.14670 


115 


18.80 + 0.02 




0.14889 


115 


18.82 + 0.02 




0.15109 


115 


18.80 + 0.02 




0.15330 


115 


18.84 + 0.03 




0.15551 


115 


18.84 + 0.02 




0.15762 


115 


18.84 + 0.02 




0.15985 


115 


18.84 + 0.02 




0.16216 


115 


18.85 + 0.02 




0.16429 


115 


18.84 + 0.02 




0.16648 


115 


18.84 + 0.02 




0.16867 


115 


18.82 + 0.02 




0.17096 


115 


18.84 + 0.02 




0.17318 


115 


18.88 + 0.02 




0.17541 


115 


18.91+0.02 




0.17765 


115 


18.95 + 0.02 




0.17998 


115 


18.97 + 0.02 




0.18214 


115 


19.01 + 0.02 




0.18437 


115 


19.03+0.02 




0.18660 


115 


19.02 + 0.02 




0.18891 


115 


19.03 + 0.02 




0.19115 


115 


19.02 + 0.02 




0.19333 


115 


19.02 + 0.02 




0.19552 


115 


19.02 + 0.02 




0.19781 


115 


19.07 + 0.02 




0.19996 


115 


19.09 + 0.02 




0.20214 


115 


19.12 + 0.02 




0.20433 


115 


19.13+0.02 




0.20667 


115 


19.15 + 0.02 




0.20886 


115 


19.15 + 0.02 




0.21105 


115 


19.16 + 0.02 





0.21324 


115 


19.18 + 0.02 




0.21556 


115 


19.20 + 0.02 




0.21767 


115 


19.21+0.02 




0.21985 


115 


19.21+0.02 




0.22204 


115 


19.18 + 0.02 




0.22433 


115 


19.21+0.02 




0.22648 


115 


19.20 + 0.02 




0.22867 


115 


19.21+0.02 




0.23086 


115 


19.21+0.02 




0.23313 


115 


19.24 + 0.03 




0.23526 


115 


19.23 + 0.02 




0.23749 


115 


19.21+0.03 




0.23977 


115 


19.21+0.03 




0.34860 


675 


19.67 + 0.02 




0.44872 


688 


19.85 + 0.02 




1.20917 


675 


21.47 + 0.06 




1.30936 


688 


21.51+0.06 




1.49485 


688 


21.64 + 0.10 




2.22673 


1727 


22.03 + 0.08 




3.15788 


1722 


22.53 + 0.13 


GROND, z 


0.12693 


115 


18.59 + 0.05 




0.12906 


115 


18.58 + 0.05 




0.13125 


115 


18.58 + 0.03 




0.13344 


115 


18.59 + 0.03 




0.13575 


115 


18.61+0.03 




0.13793 


115 


18.65 + 0.03 




0.14012 


115 


18.63 + 0.03 




0.14231 


115 


18.65 + 0.03 




0.14457 


115 


18.64 + 0.03 




0.14670 


115 


18.66 + 0.03 




0.14889 


115 


18.70 + 0.03 




0.15109 


115 


18.71+0.03 




0.15330 


115 


18.70 + 0.03 




0.15551 


115 


18.68 + 0.03 




0.15762 


115 


18.71+0.03 




0.15985 


115 


18.71+0.03 




0.16216 


115 


18.68 + 0.03 




0.16429 


115 


18.67 + 0.03 




0.16648 


115 


18.69 + 0.03 




0.16867 


115 


18.68 + 0.03 




0.17096 


115 


18.70 + 0.03 




0.17318 


115 


18.74 + 0.03 




0.17541 


115 


18.80 + 0.03 




0.17765 


115 


18.82 + 0.03 




0.17998 


115 


18.82 + 0.03 




0.18214 


115 


18.87 + 0.03 




0.18437 


115 


18.90 + 0.03 




0.18660 


115 


18.87 + 0.03 




0.18891 


115 


18.89 + 0.03 




0.19115 


115 


18.88 + 0.03 




0.19333 


115 


18.88 + 0.03 




0.19552 


115 


18.89 + 0.02 




0.19781 


115 


18.90 + 0.03 




0.19996 


115 


18.91+0.03 




0.20214 


115 


18.94 + 0.03 




0.20433 


115 


18.93 + 0.03 




0.20667 


115 


18.99 + 0.03 




0.20886 


115 


19.03 + 0.03 




0.21105 


115 


19.04 + 0.03 




0.21324 


115 


19.01 + 0.03 




0.21556 


115 


19.06 + 0.03 




0.21767 


115 


19.05 + 0.03 




0.21985 


115 


19.08 + 0.03 




0.22204 


115 


19.09 + 0.03 




0.22433 


115 


19.07 + 0.03 




0.22648 


115 


19.08 + 0.03 




0.22867 


115 


19.10 + 0.03 




0.23086 


115 


19.12 + 0.03 




0.23313 


115 


19.06 + 0.03 




0.23526 


115 


19.05 + 0.03 




0.23749 


115 


19.07 + 0.03 




0.23977 


115 


19.07 + 0.03 




0.34860 


675 


19.51+0.02 




0.44872 


688 


19.70 + 0.03 




1.20917 


675 


21.24 + 0.07 




1.30936 


688 


21.32 + 0.08 




1.49485 


688 


21.61+0.12 




2.22673 


1727 


22.00 + 0.09 




3.15788 


1722 


22.64 + 0.20 
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GROND, J 


0.13048 


730 


18.34 + 0.03 




0.13932 


730 


18.38 + 0.03 




0.14812 


730 


18.40 + 0.03 




0.15688 


730 


18.45 + 0.03 




0.16570 


730 


18.38 + 0.03 




0.17460 


730 


18.52 + 0.03 




0.18359 


730 


18.58 + 0.03 




0.19251 


730 


18.58 + 0.03 




0.20136 


730 


18.68 + 0.03 




0.21024 


735 


18.73 + 0.03 




0.21909 


727 


18.82 + 0.03 




0.22789 


730 


18.77 + 0.03 




0.23674 


730 


18.79 + 0.03 




0.34889 


730 


19.17 + 0.03 




0.44903 


730 


19.37 + 0.04 


GROND, H 


0.13048 


730.08 


18.19 + 0.05 




0.13932 


730.08 


18.18 + 0.04 




0.14812 


730.08 


18.24 + 0.04 




0.15688 


730.08 


18.27 + 0.03 




0.16570 


730.08 


18.22 + 0.04 




0.17460 


730.08 


18.31+0.03 




0.18359 


730.08 


18.34 + 0.03 




0.19251 


730.08 


18.38 + 0.03 




0.20136 


730.08 


18.51+0.03 




0.21024 


734.62 


18.53 + 0.04 




0.21909 


726.82 


18.60 + 0.04 




0.22789 


730.08 


18.55 + 0.04 




0.23674 


730.08 


18.57 + 0.03 




0.34889 


730.08 


19.04 + 0.09 




0.44903 


730.08 


19.22 + 0.09 


GROND, K 


0.13048 


730 


17.95 + 0.07 




0.13932 


730 


18.06 + 0.07 




0.14812 


730 


17.96 + 0.06 




0.15688 


730 


18.10 + 0.06 




0.16570 


730 


18.08 + 0.06 




0.17460 


730 


18.17 + 0.06 




0.18359 


730 


18.25 + 0.06 




0.19251 


730 


18.21+0.06 




0.20136 


730 


18.37 + 0.06 




0.21024 


735 


18.32 + 0.06 




0.21909 


727 


18.35 + 0.06 




0.22789 


730 


18.32 + 0.06 




0.23674 


730 


18.42 + 0.06 




0.34889 


730 


18.83 + 0.07 




0.44903 


730 


19.07 + 0.08 


F0RS2,R^;>^^^1 


0.12384 


20 


18.85 + 0.03 




0.13707 


20 


18.90 + 0.03 




0.17257 


30 


19.03 + 0.03 




0.18747 


30 


19.16 + 0.03 




0.22093 


30 


19.36 + 0.03 




0.25426 


30 


19.37 + 0.03 




0.44461 


30 


20.01 + 0.03 




1.12910 


60 


21.33 + 0.04 




1.16363 


60 


21.43 + 0.03 




1.19748 


60 


21.46 + 0.03 




1.34318 


60 


21.45 + 0.03 




1.37520 


60 


21.56 + 0.03 




1.43264 


60 


21.71+0.03 




2.36948 


120 


22.40 + 0.05 



X-Shooter, 7?^^) 



0.21572 
0.27302 



60 
60 



18.70 + 0.10 
19.01+0.09 



Gemini GMOS, r 


6.26922 


900 


23.25 + 0.04 




21.28034 


900 


23.42 + 0.05 




23.39284 


720 


23.40 + 0.10 


FT-S, Bessel R^^) 


0.66238 


1800 


20.82 + 0.02 




0.75856 


1800 


20.99 + 0.03 




0.82398 


1800 


21.23 + 0.03 




3.71735 


3600 


23.03 + 0.09 


FT-S, Sloan i 


0.63924 


1800 


20.49 + 0.02 




0.73542 


1800 


20.71+0.03 




0.80083 


1800 


20.79 + 0.05 




0.88065 


1800 


20.96 + 0.04 




3.77179 


900 


22.38 + 0.42 




3.78300 


900 


22.42 + 0.35 



UVOT white 0.0017 


149.8 


15.33+0.02 


0.0066 


19.8 


16.70 + 0.05 


0.0086 


19.7 


16.97+0.05 


0.0108 


149.7 


17.24 + 0.03 


0.0136 


19.8 


17.45 + 0.06 


0.0693 


199.8 


19.07 + 0.04 


0.1471 


361 


19.75+0.05 


0.3406 


1614.2 


20.59 + 0.07 


1.0156 


566.6 


21.73+0.18 


1.4510 


40789.7 


22.66 + 0.27 


1.9175 


29362 


22.86 + 0.19 


2.5761 


29994.8 


23.35 + 0.28 


3.2865 


46836 


22.82 + 0.18 


3.9556 


58829.1 


22.64 + 0.25 


5.2311 


138900.4 


23.79 + 0.36 


6.6682 


98684.1 


23.26 + 0.23 


8.4084 


121624.7 


> 23.81 


UVOT V 0.0072 


19.8 


16.34 + 0.14 


0.0092 


19.8 


16.84 + 0.18 


0.0132 


194.1 


16.92 + 0.14 


0.0740 


199.8 


18.51+0.17 


0.2803 


476 


19.45+0.23 


0.7414 


36383.3 


20.63 + 0.26 


UVOT b 0.0063 


19.7 


16.20 + 0.08 


0.0083 


19.7 


16.45+0.09 


0.0143 


191.8 


17.22 + 0.09 


0.1304 


1911.2 


19.38+0.14 


0.1432 


299.8 


19.19 + 0.12 


0.3259 


906.9 


20.15+0.15 


0.4128 


686.9 


20.15 + 0.17 


0.6129 


790.9 


20.27+0.18 


0.9072 


16951.2 


21.24 + 0.34 


1.0088 


603.3 


20.75 + 0.3 


UVOT u 0.0047 


249.7 


16.05 + 0.03 


0.0080 


19.7 


16.71+0.08 


0.0140 


193.5 


17.27 + 0.08 


0.0811 


199.7 


19.00 + 0.08 


0.2119 


717.9 


19.89 + 0.08 


0.4035 


907 


20.57+0.11 


0.4797 


693 


20.56 + 0.13 


0.6411 


7484.4 


21.24 + 0.14 


0.8418 


7481.9 


22.23 + 0.35 


1.5784 


63376.6 


22.69 + 0.32 


2.2041 


34276.6 


22.35+0.35 


2.9722 


87931 


22.96 + 0.3 


3.6073 


12000.2 


22.40 + 0.32 


5.0600 


132684.3 


23.01+0.35 


7.4991 


278015.6 


23.33 + 0.28 


UVOT uvwi omii 


19.8 


17.41 +0.13 


0.0098 


19.8 


17.70 + 0.15 


0.0137 


193.7 


17.83+0.11 


0.0788 


199.7 


19.80 + 0.13 


0.1885 


434.9 


20.27+0.12 


0.3930 


899.7 


20.96 + 0.13 


0.4704 


899.8 


21.35+0.17 


0.6312 


7587.9 


22.01+0.19 


0.9308 


24671.6 


22.34 + 0.24 


UVOT uvml 0.0095 


19.8 


17.66 + 0.17 


0.0124 


19.7 


18.05+0.21 


0.0764 


199.8 


20.02 + 0.17 


0.4599 


899.8 


21.29+0.17 


0.5460 


790.2 


21.70 + 0.23 


0.7027 


8389.3 


22.10 + 0.21 


0.8646 


899.7 


22.04 + 0.28 


UVOT uvw2 0.0069 


19.8 


17.88+0.15 


0.0089 


19.8 


18.05+0.16 


0.0139 


19.8 


18.52 + 0.20 


0.0717 


199.7 


20.37+0.16 


0.2581 


748.1 


21.27+0.14 


0.5256 


899.8 


21.81+0.18 


0.8314 


19056.1 


23.05 + 0.32 


1.8138 


102925.1 


> 21.90 


3.5897 


191453.2 


> 22.51 


7.1639 


325158.9 


> 22.50 



